CONSTRUCTION OF WHISKERS FOR THE QUASIPERIODICALLY 

FORCED PENDULUM 



MIKKO STENLUND 

Abstract. We study a Hamiltonian describing a pendulum coupled with several 
anisochronous oscillators, giving a simple construction of unstable KAM tori and their 
stable and unstable manifolds for analytic perturbations. 

We extend analytically the solutions of the equations of motion, order by order in 
the perturbation parameter, to a uniform neighbourhood of the time axis. 



1. Main Concepts and Results 

1.1. Background and history. A quasiperiodic motion of a mechanical system is 
composed of incommensurable periodic motions; the trajectory in phase space winds 
around on a torus filling its surface densely. An integrable Hamiltonian system has a 
great profusion of quasiperiodic motions: if one picks an initial phase point according to 
a uniform distribution, the trajectory will be quasiperiodic with probability one. The 
remaining trajectories are periodic. 

KAM theory deals with the stability of quasiperiodic motions, or persistence of invari- 
ant tori, under small perturbations. Poincare [Poi93a] called this the general problem 
of dynamics. 

In 1954, Kolmogorov [Kol54] outlined a result, made rigorous by Arnold in 1963 
[Arn63], that quasiperiodic motions are typical also for nearly integrable analytic Hamil- 
tonians under suitable nondegeneracy conditions. Thus, only a small fraction of the tori 
would be destroyed by the perturbation. Moser managed to prove the same for twist 
maps [Mos62] in 1962, and later for Hamiltonians [Mos66b, Mos66a], in the smooth 
(non-analytic) setting (see also [Mos67]). 

The difficult problem to overcome is the following. Suppose that the Hamiltonian 
reads 7i = Tio + XHi, where H,q is integrable and A is considered small. Then one can 
formally represent a solution to the equations of motion by a power series in A, known 
as the Lindstedt series in this context, conditioned to agree for A = with a quasiperi- 
odic solution obtained in the integrable case. When one computes the coefficients of 
the Lindstedt series, however, one encounters expressions containing arbitrarily small 
denominators. The latter seem to imply that the fcth coefficient grows like k\ a with a 
large power a. Thus, there is little hope of being able to sum the series and obtain a 
true solution, unless a miracle occurs. 

The proofs mentioned above relied on a rapidly convergent Newton-type iteration 
scheme, which is interesting in its own right, and yields solutions analytic in A. On the 
other hand, one is then left to wonder why the Lindstedt series does converge. 
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In 1988, an answer was provided by Eliasson [EH96], who managed to identify enor- 
mous cancellations among the small denominator contributions and to sum the Lind- 
stedt series "manually" . Gallavotti [Gal94a, Gal94b] interpreted the cancellations in a 
Renormalization Group (RG) framework. For a review and some extensions, see Gentile 
and Mastropietro [GM96]. The importance of these achievements has to be stressed: 
they prove the existence of quasiperiodic solutions in an essentially constructive way. 

Motivated by the RG approach of Gallavotti, in the 1999 paper [BGK99] Bric- 
mont, Gaw§dzki, and Kupiainen identified the cancellations as a consequence of Ward 
identities (corresponding to a translation invariance of an action functional) in a suitable 
field theory. 

Returning to much earlier works, Moser [Mos67] and Graff [Gra74] showed that also 
hyperbolic tori — tori having local stable and unstable manifolds — would typically persist 
under small perturbations. In another landmark paper [Arn64] , Arnold had described a 
mechanism how a chain of such "whiskered" tori could provide a way of escape for special 
trajectories, resulting in instability in the system. (As discussed above, a trajectory 
would typically lie on a torus and therefore stay eternally within a bounded region 
in phase space.) The latter is often called Arnold mechanism and the general idea of 
instability goes by the name Arnold diffusion. It is conjectured in [AA68] that Arnold 
diffusion due to Arnold mechanism is present quite generically, among others in the 
three body problem. 

Arnold mechanism is based on Poincare's concept of biasymptotic solutions, discussed 
in the last chapter of [Poi93b], that are formed at intersections of whiskers of tori. Fol- 
lowing such intersections a trajectory can "diffuse" in a finite time from a neighbourhood 
of one torus to a neighbourhood of another, and so on. 

Chirikov's work [Chi79] is a very nice physical account on Arnold diffusion. Lochak's 
compendium [Loc99] discusses more recent developments in a readable fashion and is a 
good point to start learning about diffusion. 

The proofs of Moser and Graff mentioned above use the rapidly convergent method 
of Kolmogorov, but there now exist also constructive proofs in the spirit of Eliasson and 
Gallavotti. We refer here to Gallavotti [Gal94b] and Gentile [Gen95a, Gen95b]. 

1.2. The model. We consider the Hamiltonian 



of a pendulum coupled to d rotators, with G S 1 := R/27rZ and / el the coordinate 
and momentum of the pendulum, and ip G T d := (S x ) d and A G M. d the angles and 
actions of the rotators, respectively. The perturbation / is assumed to be real- valued and 
real-analytic in its arguments, and A is a (small) real number, whereas the gravitational 
coupling constant g is taken to be positive. This Hamiltonian is sometimes called the 
generalized Arnold model or the Thirring model. It is the prototype of a nearly integrable 
Hamiltonian system close to a simple resonance, as is explained in the introduction 
of [Gen95b]. A review of applications can be found in [Chi 79]. 
The equations of motion are 



H(1, 0, A, = \l 2 + g 2 cos + \A 2 



(1.1) 



= /, ijj = A, I = g 2 sm</) + Xd^f, A = Xd^f. 



(1.2) 
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For the parameter value A = 0, which is addressed as the unperturbed case, the pen- 
dulum and the rotators decouple. The former then has the separatrix flow : R — > S 1 
given by 

<P{t) = $ V), 

where 

$°(z) = 4arctanz. 

By elementary trigonometry, this function possesses the symmetry property 

$°(z) = 27T-$°(^ 1 ). (1.3) 

It is also odd, 

<f>°(-z) = -$°(z). 

The phase space of the unperturbed pendulum looks as in Figured! where the separatrix, 
given by <3>°, separates closed trajectories (libration) from open ones (rotation). 




Figure 1 . A (0, J) plot showing the unperturbed pendulum separatrix 
that intersects the axis at integer multiples of 2n — the upright position 
of the pendulum. 

On the other hand, ip : R —>■ T d is quasiperiodic: 

il>(t)=ij)(p)+ut (mod 2tt), 

such that the vector 

u := A(0) = A(t) 

satisfies the Diophantine condition 

\oj ■ q\ > a\q\~ v for q G Z d , q ^ 0, (1.4) 

with a and v positive. Thus, at the instability point of the pendulum, the flow possesses 
the invariant tori 



{(0,^/,A) = (o,M,ccO|#eT d } 



indexed by lu, with stable and unstable manifolds (Wq and Wq , respectively) coinciding: 

Wr = {(M,I,A) = (<£ (z),6,gzd z $°(z),uj) | z E [-00,00], BET*}. (1.5) 

Remark 1.1. The constant g is the Lyapunov exponent for the unstable fixed point of 
the pendulum motion; in the limit s — > —00 two nearby initial angles <j)(s) and <f)(s + 5s) 
separate at the exponential rate e gs . As <p(t) = $°(e t//ff ), the Lyapunov exponent 
fixes a natural time scale of g~ x units, characteristic of the pendulum motion in the 
unperturbed Hamiltonian system (II. ip . 
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When the perturbation is switched on (A 7^ 0), we show that some of the invariant 
tori survive and have stable and unstable manifolds — or "whiskers" as Arnold has called 
them — that may not coincide anymore. 

1.3. Main theorems. Our approach will be to construct the perturbed manifolds in a 
form similar to (jl.5p as graphs of analytic functions over a piece of [—00, 00] x T d . To 
see how this can be achieved, note that the unperturbed stable and unstable manifolds, 
Wq and Wq , consist of trajectories 



that at time ±00 become quasiperiodic, as they wrap tighter and tighter around the 
invariant torus %; indeed (4>(t),ip(t)) ~ (0, cut) in the limit t — > ±00. 

Analogously, we will find the stable and unstable manifolds of the perturbed tori by 
looking for solutions of the form 



with quasiperiodic behavior in one of the two limits t — > ±00. Note especially that we 
anticipate the Lyapunov exponent 7 > to depend on A, with 7| A=0 = g. 

Remark 1.2. One should not assume asymptotic quasiperiodicity in both of the limits 
t — > ±00, as the unstable and stable manifolds, which we denote W" and W|, are 
generically expected to depart for nonzero values of the perturbation parameter A. 
Therefore, either the past or future asymptotic of a trajectory will evolve so as to 
ultimately reach the (deformed) invariant torus T\. The separatrix in Figured] is thus 
transformed into something like the pair of curves in Figure [2J 




Figure 2. A schematic i-versus-0 plot, on a section of constant ip (d = 
1). The stable and unstable manifolds are expected to split, as opposed 
to coincide. The origin has been shifted for convenience. 

Let us denote the total derivative d/dt by d t and the complete angular gradient 
(d,!,, by d for short. Substituting (11.61) into the equations of motion 




(<i>(t),i/>(t)) = ($(e 7i ,^),^ + ^(e 7i ,wt)) = (0,c^) + ($,^)(e 7t ,c^) 



(1.6) 
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Notice that the partial differential operator 

£ : = uj ■ d e + ^zd z 

satisfies the characteristic identity 

CF(ze yt , 9 + ut) = d t F(ze lt , 9 + ut) (1.7) 

for a differentiable map (z, 9) i— > F(z,6). Equation (11.71) simply reflects the time deriv- 
ative nature of C In fact, if T is the "time-reversal map" 

T( z ,e) = (z-\-e), (i.8) 

then, by the chain rule, 

C(FoT) = -(CF) oT. (1.9) 

Let us abbreviate 

Q(X) :=(g 2 sm^,0) + XQ(X) with Q(X) := df(X + (0,9)). (1.10) 

As a consequence, we have reduced the equations of motion to the PDE 

C 2 X = Q(X) (1.11) 

for the map (z, 9) X(z, 9) in a suitable Banach space of analytic functions, albeit its 
restriction to the set ("characteristic") 

{(z,9) = (e~<\ujt) | t G R} (1.12) 

is what one is physically interested in. Our preference of working directly with the in- 
variant manifolds, as opposed to individual trajectories traversing along them, motivates 
us encoding the time derivative in the operator C Nevertheless, it will be harmless — 
and indeed quite informative — for the reader to keep in mind that the objects we deal 
with originate from (11.121) and therefore have a direct physical interpretation. 
The action variables trivially follow from the knowledge of X(z, 9): 

(I(t),A(t)) = (0,u) + Y(e rt ,ut), Y := CX. 

The solutions X will provide a parametrization of the deformed tori and their stable 
and unstable manifolds. As hinted below (II. 6p . we find two kinds of solutions, X u (z, 9) 
defined for z G [— z , z ] ='■ I" and X s (z, 9) defined for z G [— oo, — z^ 1 } U [zq 1 , oo] =: P. 
Here, zq > 1. The tori will have the three parametrizations 

r A = { (^,i,A) = ((o,0) + x u (o,e), (o,w) + 1^(0,0)) I eeT d ] 

= {(^^,I,A) = ((0,9) + X s (±^ i e) 1 (0,u) + Y s (±^,9)) \ 9eT d }, 
whereas the parametrizations of their stable and unstable manifolds then read 

W s x u = { (</>, V, /, A) = ((0, ff) + X s ' u (z, 6), (0, w) + Y s - U (z, 6))\ze I s ' u , 6eT d }. 

In order to enable solving (11.111) . we need to deal with quantities of the form (uo ■ g) _1 , 
q G Z d \ {0}, stemming from the Fourier representation of the operator C. Here the 
Diophantine property of the vector uj G M. d stated in (11.41) steps in. Since u = A\ x=0 = 
-0| A=O , by rescaling time (and the actions, correspondingly) in the equations of motion 
(II .2p . the constant a can be absorbed into g 2 and A in the equations of motion, leaving 
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the ratio \g~ 2 unchanged: (g, A) i— > (g/a, A/a 2 ) Q. Thus, we may as well take a to be 1 
below, transforming the condition on uj into 

\uj-q\ > \q\~ v for qeZ d \{0}. (1.13) 

We will moreover consider A small in a ^-dependent fashion, taking 

e:=\g- 2 (1.14) 

small. This should be seen as an outreach towards the experimenter, albeit there is 
a technical wherefore: such a choice is needed for studying the limit g — > oo, which 
corresponds to rapid forcing; see Remark 11.11 The domain we restrict ourselves to is 

D := {{e,g) E C x R \ \e\ < e , < g < g } , (1.15) 

for some positive values of e and g . 

Finally, note that if X = ($, \P) solves (11.111) on some domain D' C [— oo, oo] x T d , 
then so does 

X a> p(z, 9) := X(az, 9 + 0) + (0, /3), (1.16) 

as long as (az, 9 + (3 (mod 2ir)) G D'. The aforementioned invariance is a manifestation 
of the freedom of choosing initial conditions for (0, ip) — we may choose the origin of 
time and the configuration of the physical system there. 
For e = 0, the solutions are obtained from 

X°M):=($°(z),0) (1.17) 

using (I1.16p . In particular, X°(1,0) = (tt,0). This will provide us with a natural way 
of fixing a and j3 below. 

We are now ready to state the first of the two main theorems of this article. It is 
a version of a classical result, and by no means new; earlier treatments include for 
instance [Mel63, Mos67, Gra74, Eli94, Gal94b, Gen95a, Gen95b]. However, the interest 
here lies in the new techniques used in the proof. 

Theorem 1 (Tori and their whiskers). Let f be real-analytic and even, i.e., 

Also, suppose uj satisfies the Diophantine condition (11.131) . and fix g$ > 0. Then there 
exist a positive number 6q and a function 7(e,g) on D, analytic in e with \j — g\ < Cg\e\, 
such that equation (11. lip has a solution X u which is analytic in e as well as in (z, 9) in 
a neighbourhood of [—1, 1] x T d and which satisfies 

X"(1,0) = (vr,0), X u (z,9) =X (z) + O(e). (1.18) 

Corresponding to the same 7, there exists a solution X s (z,9) = X°(z) + 0(e) which is 
an analytic function of {z^ 1 , —9) in a neighbourhood of [—1,1] x T d . The maps 

W s ' u {z,9) = {X s ' u ,Y s ' u ){z,9) + {{0,9),{0,u)), Y s ' u := CX S ' U , (1.19) 

provide analytic parametrizations of the stable and unstable manifolds W^'" of the torus 
T x . 

^This scaling is responsible for the usual requirement A = 0(a 2 ) for KAM tori. 
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Remark 1.3. The number eo above depends on the Diophantine exponent v and on /. 
The perturbation (4>,ip) i— > f(<p,ip) is analytic on the compact set S 1 x T d . By Abel's 
Lemma (multivariate power series converge on polydisks), it extends to an analytic map 
on a "strip" |3tn</>|, |9m-?/>| < r\ [r\ > 0) around S 1 x T d . By Theorem [fl there exists 
some < a < i] such that each 9 t— > X s,u ( ■ ,6) is analytic on |9m#| < a. 

An important part of Theorem[T]is that the domains of X u and X s overlap. Namely, if 
0,0) i — ► X(z,6) solves equation (jTTTl) . then so does (z,6) h-> (2ir,0)-(XoT)(z,9). This 
is due to (11. 9p and the parity of /. Consequently, by a simple time-reversal consideration 
(set t i — ► — t in ( I1.12p ). the stable and unstable manifolds are related through 

X s = (27r, 0) — X u o T. (1.20) 

In particular, as T(1,0) = (1,0), 

A S (1,0) = A"(1,0). 
Moreover, the actions Y s,u = CX S,U satisfy 

Y s = Y u oT 1 (1.21) 

yielding 

F S (1,0) = F M (1,0). 

In other words, a homoclinic intersection of the stable and the unstable manifolds W^' M 
occurs at (z, 9) = (1,0), as their parametrizations (11.191) coincide at this homoclinic 
point. Since the manifolds W^' n are invariant, there in fact exists a homoclinic trajectory 
on which the parametrizations agree: 

iy s (e 7i , out) = W u {e^\ tot). (1.22) 

Remark 1.4. Equation (11.201) is what remains of the symmetry X° = (2tt, 0) — X° o T, 
which is just another way of writing fll.3p . after the onset of even perturbation. This is 
an instance of spontaneous symmetry breaking: The equations of motion, (II. lip , remain 
unchanged under the transformation X i— > (2tt, 0) — X o T, but the individual solutions 
do not respect this symmetry; X u ^ X s = (2ir, 0) — X u o T, if A ^ 0. 

Coming to the second one of our main results, let us expand 

oo 

x u = j2 ^ xu/ - 

In Section we will show that the common analyticity domain of each X u ' in the 
z- variable is in fact much larger than the (small) neighbourhood of [—1, 1] — the cor- 
responding analyticity domain of X u according to Theorem [lj namely it includes the 
wedgelike region 

U Ti1? := {|*| <r} |J{arg2G [-0, 0] U [tt - tx + 0] } CC 

(with some positive r and 0): 
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Theorem 2 (Analytic continuation). Each order X u ' of the solution extends ana- 
lytically to a common region x {|$5m#| < a}. Moreover, if if> i— > f(- ) %p) is a 
trigonometric polynomial of degree N , i.e., N is the minimal nonnegative integer such 
that f(-,q) = whenever \q\ > N, then 9 \— > X u ' e ( ■ , 9) is a trigonometric polynomial 
of degree £N, at most. 

Remark 1.5. With rj and a as in Remark 11.31 the numbers r and •& are specified by 
the following observation: $°(V) = 4arctanz implies that |Q ; m < l )CI (z)| < r] in U T! i? 
with r and $ sufficiently small. By Remark 11.31 ( z i@) l— ^ f{^°{ z ),@) is analytic on 
U T) # x {l^m^l < a}, which we will use as the basis of the proof. 

In spite of Theorem EJ (a straightforward upper bound on) X u ' grows without a limit 
as |9fte z\ — > oo, such that there is no reason whatsoever to expect absolute convergence of 
the series YH>Lo £ £ X u,e in an unbounded z-domain with a fixed e. In fact, it is known that 
the behavior of the unstable manifold gets extremely complicated for large values of z 
even with innocent looking Hamiltonian systems. Still, it seems to us that the possibility 
of a uniform analytic extension of the coefficients X u,e has not been appreciated in the 
literature. 

Due to fll.20p , an analog of Theorem [2] and the subsequent discussion are seen to hold 
for the solution X s , with z replaced by z~ l . 

Theorem[2]is interesting, because it allows one (at each order in e) to track trajectories 
t I— > W s ' u (e' yt , 9 + ut) on the invariant manifolds W^'" for arbitrarily long times in a 
uniform complex neighbourhood |Qtn.t| < g~ lq & of the real line, for arbitrary 9 G T d . 
The motivation for doing this stems from studying the splitting of the manifolds W^' M 
in the vicinity of the homoclinic trajectory (11.22p . and is the topic of another article. 
The general ideology that, being able to extend "splitting related functions" to a large 
complex domain yields good estimates, is due to Lazutkin [Laz03], as is emphasized 
in [LMS03] . 

1.4. Strategy. Let us briefly explain how Theorem [1] will be proved in three steps. Due 
to fll.20p . we may concentrate on studying the unstable manifold. Thus, we write 

X(z, 9) := X u (z, 9) = X {6) + zXxiO) + S 2 X(z, 9). 

From (11. lip we first get an equation for Xo := X u (0, •) alone. Second, given Xq, an 
equation for X% := d z X u (0, ■ ) and 7 alone is obtained. Third, given X , X\, and 7, an 
equation for the remainder 8 2 X is obtained. 

It turns out that solving for X and X% (together with 7), i.e., the invariant torus and 
the linearization of the unstable manifold around it, is difficult. Namely, these problems 
involve the small denominators of KAM theory. In contrast, solving for 8 2 X amounts 
to a simple Contraction Mapping argument. 

We deduce the existence of Xq from [BGK99] . The existence proof of X\ is reminiscent 
of the RG argument in the latter paper, except that the Lyapunov exponent 7 has to 
be fine-tuned to a proper value such that the renormalization flow converges. 

At this point we would like to draw the readers attention to the interesting refer- 
ence [Gen95b] , where the author takes a different approach. Gentile fixes the perturbed 
Lyapunov exponent 7 in advance and replaces g by fif(e, 7) in the Hamiltonian, which 
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is analogous to introducing counterterms in quantum field theory, and finds the cor- 
responding manifolds. One could then solve the implicit equation g(e,j) = g and to 
obtain 7 as a function of g and e. 

Acknowledgements. I am indebted to Antti Kupiainen for his help during the course 
of this work. Guido Gentile, Kari Astala, and Jean Bricmont provided sharp remarks 
and critical comments on the manuscript that made it more comprehensible and math- 
ematically accurate. I wish to express my gratitude to all of them. I thank Giovanni 
Gallavotti and Emiliano De Simone for discussions at Rutgers University and University 
of Helsinki, respectively. 

2. Perturbed Tori 
The perturbed tori will be found by looking for solutions having the general form 

<j)(t) = %(uot), ^(t)=ut + %(ut), 
with $ : T d -> R and ^ :T d ^ R d satisfying the "t -»• -00 asymptotics" 

V 2 <!> (6) = g 2 sin $ o (0) + A fy/($ o (0), + %(9)) (2.1) 
V 2 V o (6) = \d^f(%(6),6 + M0)) (2-2) 

obtained from equation (11.111) by putting z = and V = to ■ d$. Note that if X = 
($0, ^0) is a solution to equations (12.11) and (12. 2p . then so is 

apXoiff) := + P), *o(0 + 13) +13) (2.3) 

for (3 E T d . We point out that together (12.11) and (12.21) are equivalent to 

V 2 X = n(X ). (2.4) 

2.1. Spaces of analytic functions. Let us define the spaces we shall be working in. 
As linear subspaces of i 1 , the Banach spaces 

Bf := |$ : T d -> C := \Ho) < °°}' 

q&Z d 

Bl := I* : Jd ^ C d := ^|^(g)|e ff|g| < 00}, 



for any a > 0, have the advantage that Fourier analysis on their elements is convenient. 
Furthermore, we are trying to find a solution X = ($, \l/) analytic on the torus, and, for 
a suitably small a, such a function belongs to B® x £>* because of the exponential decay 
of its Fourier coefficients; |X(g)| < Ce~°"' 9 ' with some positive constant C. Indeed, if 
a > 0, the spaces above comprise precisely those functions on the torus that admit 
an analytic extension to the "strip" |9tot0| < o. We will occasionally write B a when 
referring to either one of B® and . 

Of course, as our analysis proceeds, the perturbation / will appear all over the place. 
This, in turn, dictates the analyticity properties of a plethora of maps, in practice 
introducing the constraint a < 77 for the spaces B a ; see Remark 11.31 
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Notice the natural embeddings 
for a > 0, due to the inequality 



B a +a C B a , 



< II ■ IL„- (2-5) 



Consider the linear operator Tg : B a+a — > B a defined through setting TpX(q) = 
e n-f>X{q), with (3 G C d . Whenever |9fm/3| < a, IMI^ Bff) < 1. The realization of 
Tp in terms of the variable 9 is just the translation ^(6) i— > ^(0-1-/3). will serve as 
a useful device in encoding the real-analyticity of / as an algebraic property into the 
Fourier series of certain other functions. This is due to the the fact that exponential 
smallness of |X(g)| in q implies real-analyticity of a function X on the torus, and vice 
versa. 

We shall encounter n-linear maps from C d+1 into C. Endowed with the norm 

H A ll£("(c*+i);C) :=mf{M>o| \A{z u . . . , z n )\ < M\ Zl \ . . . \z n \ Vz t e C^ 1 } 
they form the Banach space £( n (C d+1 ); C); see [Cha85]. 

2.2. Past and future asymptotics of the solution in the perturbed case. This 
subsection discusses the t — > ±oo asymptotics of the solution X. In these limits the mo- 
tion settles onto the "distorted version" T\ of the invariant torus T with the pendulum 
seizing to swing, but wiggling quasiperiodically about its unstable equilibrium. 

Theorem 2.1. Under the assumptions of Theorem^ there exist positive numbers r 
and €q such that, for (e, g) G D, equations (12.1 1) and (12. 2p have a unique solution X = 
($0,^0) i n the class of those real- analytic functions of 6 G T d that satisfy H^olU < t 
and (^o) = (zero average). The function Xq, defined on {\Qm8\ < a} x D for some 
a > 0, is analytic and uniformly bounded by (C\e\, Cg 2 \e\). Moreover, it is ~RxM. d -valued 
on T d for e real. Thus, any real-analytic solution X' Q = (& ,ty ) with (ty' ) = (3 G R d 
and \\^' — /3\\ £l < r must be the one given by 

x' (e) = x (e + ^ + (o^), 

i.e., Xq = o~pX , using the notation of (12.31) . 

Remark 2.2. Remark 11.31 below Theorem [TJ holds true. Recall that we have defined 
e := Xg~ 2 in fll . 14j) and the domain D in (I1.15p . This is a version of the KAM Theorem. 
Notice that X G B* x B* . 

Proof. The proof is a reduction to the one given in [BGK99]. Here we systematically 
omit the subindex of $ , an d X . Let us concentrate on the pendulum part, 
equation (12.11) . first. We expect $ to be close to its unperturbed value, zero, and it pays 
to cancel the leading term of g 2 sin$(#) on the right-hand side by subtracting g 2 §(6) 
from both sides. We then have 

(V 2 -g 2 )<$> = U(<S>,y) =: U(X) (2.6) 

with 

U{X)(B) := /(sin $(0) - $(9)) + X fy/($(0), B + V(6)). (2.7) 
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Pay attention to the fact that U(X)(8) depends locally on X — only through X{9), 
that is. Abusing notation, we shall use U(X)(8), U(X,8), U{X{9),9), etc., in the 
same meaning, whichever is the most convenient form. Now, U(x, 0) is analytic in the 
vector argument x — (Xfr Xif>) m ^ ne region \x<t>\, \Xi>\ — where r\ > depends on the 
analyticity domain of /; see Remark 11.31 on page CD 

Let us now write down the Fourier-Taylor expansion 

00 1 

U(X(B), 8) = J2 ~l Dn U(0, 9) (X(9), X{9)) 

n=0 

00 

= E e i0 ^D n U(O,9)(X( qi ),...,X(q n )), (2.8) 

n=0 ' q=(gi,...,g n ) 

qi& d 

where D n U(0,9) G £( n (C d+1 );C) is the nth Frechet derivative of the map U( ■ , 9) : 
C*+i - C : X » U( X ,9). 

The map 9 \— > U n {9) := ^D n U(0, 9) is analytic in the same domain as 9 \— > U(0, 9) = 
A 9^/(0, 9), i.e., \Qm9\ < r\. Its Fourier representation U n {9) = J2 q ez d e%q 9u n(q) has 
coefficients 

= 7^T- d I e~ iq (> -D n U(0, 9) d9 (2.9) 
in £( n (C d+1 ); C). Using this notation, we translate (I2.8P into the Fourier language; 

oo n 
^X<?) = E E Wn(g-E^(^)>---^(?n)). (2.10) 

The right-hand side of equation (12.101) is a power series in X, converging whenever 
X is sufficiently close to zero. Namely, we have 

Lemma 2.3. The multilinear maps u n (q) obey the bound 

IK^il^c^c) < Cg\rl + |e|)(r /e)-"e-^, (2.11) 

where p and r is any pair of positive numbers satisfying p + r = r],r]>0 being the 
width of the analyticity domain of f as explained in Remark ] 1.31 

The proof of Lemma |2~31 is straightforward, but, for the sake of continuity, is given in 
Subsection 12.31 below. 

Considering the closed origin-centered balls of radius r < r /2 in £>* and £>* — Bf r and 
B% T , respectively— we next study Up : B* r x B% r -> : ($,*) h-> t p U{t. p ^,t^). 
By equation (12.71) . 

^($(0), M>(6>), 0) = U($(6), 8 + P + (2-12) 

when /? G M d . The right-hand side is analytic in /?, and extends to |3fm/3| + a + r < 77 
through the same expression, leaving Up analytic with respect to X. 
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More quantitatively, one checks using the bound ( 12. lip that the power series 

oo n 

W0(?)=£ E e i ^-^K n (q-J2<li)(X(<lil---iX(<ln)), (2-13) 

n=0 qe(Z d ) n i=0 

converges uniformly with respect to X and /3, even if the latter has a small imaginary 
part. In fact, || ct obeys the upper bound 

oo n 

E E E e<Tkl | e ^ E<9i) ^^-E*)(^(ft)>---»^(?n)) 

n=0qe(Z d )" ggZ d «=0 

oo / \ n 

< E E e (|3m/3|+CT)kl ii^(g)ii £( . (cd+1); c) E ni 1 ^)i eff|gi1 

n=0 \q£l d J qe(Z d ) n *=l 

oo 

< Cg\rl + kl) E E e^l+^'M (r /e)-" \\X\\ n a < Cg 2 (r 3 + |e|), 

if we choose |Qtn/?| + a < p = r\ — r and r < r /2e, since ||X|| CT < 2r. Thus, fixing 
r = r /6, say, we obtain 

sup ||^(X)|| CT <C/(r 3 + |e|) (2.14) 



whenever 

\Qmp\ + a + 6r <r]. (2.15) 

Lemma 2.4. Suppose ( I2.15P holds, and \I> G B^ r . Then, for r and eo small enough, 

(V 2 -g 2 )$ = Up($,V) 

has a solution Qp^) G Bf r , real-valued provided (3, e, and \l/ are, and there are no 
other solutions in the t x -ball B^ r d . In fact, $g(\I/) = rg$ (r_ i g\E'). The map 
^ i — > $^( X I / ) is analytic on B^ r . $ ( a(\I / ) a/so depends analytically on (3 as well as on 
(e, g) G -D fsee (11.151) ). and o6e?/s iae bound 

||^)L<C|e| (2.16) 

uniformly in ^ , {3, and g. 

Remark 2.5. The smallness condition is C(r 3 + eo) < r, where C is the same constant 
as in (I2.14p and contains the norm of the perturbation /. 

The standard but lengthy proof of Lemma 12.41 may be found in Subsection 12.31 
Let us come back to equation fl 2 . 2 j) . whose right-hand side may now be written solely 
in terms of ^ G -B* r , amounting to 

V 2 ^ = V(V) (2.17) 

with V(V)(6) = Afy/($(tf)(0),0 + V(0)). Consider then V^(*) := t p V{t_^). By 
Lemma [2 A\ it reads 

Vf,(*)(e) = V(t^)(6 + (3) = A^/(^(*)(0), 6 + (3 + *(6/)) 
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and is analytic in the domain 

B* xflx {\%m9\ <a}x{f3\ \$SmP\ + ct + 6r < 77} (2.18) 
with the uniform bound 

||W0L< sup |A^/(0,V)i<C^ 2 |e|, 

\^m<f>\,\^mip\<r) 

provided C\e\ < r\ (see (I2.16P ). 

Equation (12.171) is the variational equation corresponding to the action functional 



S : B* r -> R : * i-> = / s(tt, 0) d# 



given by the integrand 

S (*, 0) = ±($£> 2 $ + * . £>2^) + ^2 cog $ _ A y^ ; + 

where $ = 5 is invariant under the T d - action *(0) i-> ^(0) := \l/(0 + (3) + (3, 

(3 G R d . Hence, dpS(^p)\p =0 = yields the WW identity 

jCW*=jC ,(s) -*'S* (i=1 '-" d) (2 - 19) 

of the symmetry in the functional derivative notation. In fact, 

(£> 2 ^-y(#))(0). 



Integrating by parts three times one sees that 

/ V(9) ■ d e ,.V 2 ^(9) dd = - I m(d)-d e iV 2 m(6)dd = Q. 
./• •••' ./ •• 

The general identity (12.191) therefore reduces to the identity 

V i (V,6)d6= [ ^(9)-d 6 iV(^,d)de (2.20) 



jd jjd 
for the map V. 

In conclusion, we have the KAM-type small denominator problem (I2.17P with Vp{^, 9) 
analytic in the domain (12.181) and bounded there by C|A|, together with the Ward 
identity (I2.20p stemming from a translation symmetry of the action that generates 
the equation. Furthermore, Va(\I/,0) is real-valued whenever (3, e, and \I/ are. For 
< a < i] — 6r — so that we may choose ^ — this is precisely the setup in 

[BGK99], where the authors devise a method for dealing with such problems using a 
Renormalization approach. 

The subtle analysis in [BGK99] yields a unique solution ^ G B^ r to (12.171) with zero 
average and analytic in (e, g) G D. The inevitable loss of analyticity takes place in 
the domain of (3. The map 9 i— ► ^(9) is R -valued on the torus for real e and satisfies 
ll*IL<C|A| = C7^|e|. 



14 



MIKKO STENLUND 



Denote by \I/ n , n G Z+, the unique solution to ( 12.171) in the ball B^ nr . Since _B* r C 
Bp/ nr , \P has to coincide with \l/ n . Hence, \1/ is the unique solution in 

oo 

[J Bf /nr D |^ : T d -»• M d ^ real-analytic and ||*||^ < r| . 

n=l 

Indeed, assuming the map i— » \l/(#) is real-analytic, || 1 P|| a ./„ < oo for some n, and we 
have that H^H^ \ ll^llo = II^IU as n ^ oo. Thus, if H^l^i < r, we gather that 
ll^llcr/n < r ^ or sufficiently large values of n. 
This concludes the proof of Theorem 12.11 □ 

2.3. Proofs of Lemmata E3] and E31 

Proof of Lemma \2.3[ Write || ■ || = || ■ \\c( n (c d+1 yc) ^ or short. From (I2.9p and the Cauchy 
Integral Theorem, 

K(<?)H= 7^ / e-^< 9+ ^ -D n U{% 9 + i(3) d6 

[ZTT) Jjd n\ 

< e q - p — sup||L> n [/(O,0 + z/3)||, 

for f3 G W 1 and \(3\ < r\. Take < p < r\ and choose (3 = —pq/\q\- We compute the 
standard norm of n-homogeneous polynomials, 

\\D n U{0, 9 + || P „ (cd+ i. c) := sup |,D n t/(0, + 0) (z, . . . , z)\, 

\z\<l 

which, using the Cauchy Integral Formula, turns into 
n! r U((z,e + i/3)d( 

'9B(0,r ) 



sup 

\z\<l 



7T-.f Vn+i <n!r Q ~" sup sup |f/(C^ + i 

7anjfo,roi <•> |2|<iceaB(o,ro) 



Here O(0, r ) is the origin-centered circle of radius r in the complex plane, with the 
constraint r + p < rj. For \z\ < r and |3m#| < p we estimate 

<C(? 2 (r 3 + |6|); 

see equation (12. 7\i . Here we have singled out \g~ 2 = e, and C is independent of g. 
We stress that U(z, 9) simply stands for the expression obtained from the expression of 
U{X, 9) in (|2Z3> by replacing X(9) by z G C d+1 . 

Symmetric multilinear maps are fully determined by their diagonal — the correspond- 
ing homogeneous polynomial, that is — which is explicitly confirmed by the Polarization 
Formula [Cha85, Din99]. Hence, in order to obtain the estimate in ( 12.111) . we multiply 
the corresponding polynomial estimate by the factor n n /n\ ~ e n /y2nn. □ 



Proof of Lemma \2.4\ The proof is a simple application of the Banach Fixed Point The- 
orem. We fix * G B* r and study the operator F($) := (V 2 - g 2 )' 1 !/^, *). 
First, {T> 2 — g 2 )' 1 is a linear operator bounded in norm by g~ 2 . From (12.141) . 

\\F(*)\l < g- 2 \\U p (<S>,ni < C(r 3 + \e\) < r 
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for sufficiently small r and e, which means that F(B® r ) C Bf r . Proving contractiveness 
resembles estimating the norm of Up in the proof of Theorem 12.11 and is omitted. The 
existence and uniqueness of the solution <3>(\I/,/3) G B® r now follow. 

For (3, e, and \P real, F maps the closed subset of real-valued functions $ G Bf r into 
itself and is a contraction there, so $( 1 I r ,/3) is real-valued by uniqueness. 

The operator F depends analytically on the parameter \I/ in B^ r . Consider the 
sequence (F k (0)) keN of successive substitutions. Each element F k (0) is analytic in 
\P G -B* r . Furthermore, the Banach Fixed Point Theorem guarantees that such a 
sequence converges to the fixed point $( 1 I r ,/3) in geometric progression; 

||F fc (0) - *(%0)h < W F (°)L < rZT- 

Consequently, $(^ r ,/3) is the uniform limit of a sequence of analytic functions, and, as 
such, analytic itself. The same argument goes for (e,g) G D (see ()1.15p ). as well as for 
(3 in the domain specified by (12.151) . 

Because ( 12. 5p implies ^ G -B* r C -B* r , in fact $( 1 i r ,/3) is the unique solution in £>* r . 

Let us denote = $(tf,0). If * G £?* r and |3?m/3| < a/2, then r./j* G 

S*, 2r , such that $ = ^(t^^) is the unique element in B®/ 2r solving $ = (V 2 — 
g 2 )~ 1 U(&,T_ /3 ^). The diagonality of rp and £> yields 

where &p(ty) = TpQ(r_pty) G £>* r . But was the unique solution in £?* r , such 

that &p(ty) = $(^, 0) G -B* r . For larger |Sm/5| one obtains an analytic continuation. 

A priori, we know that ||<l ) ( g(\[ / ) || CT = ||-^ 1 (^ y s( 1 S r ))|| (T < t. On the other hand, we know 
that $^(^ r )| e=0 = by uniqueness, whence the estimate (12.161) follows. □ 



3. Lyapunov Exponent — Linearizing the Unstable Manifold 

In this section we study the motion in the immediate vicinity of the torus T\ cor- 
responding to the solution X (9) of Theorem 12. 1[ To that end, suppose X(z,9) is an 
analytic solution to equation (fTTTj) with X(0,6) = X (9). Then Xi(9) := d z X(0, 9) 
should satisfy the equation 

(V + 1 ) 2 X 1 =Dn(X )X u (3.1) 

as Q(X)(z, 6) depends on z only through X evaluated at (z, 6). 

Note that (13.11) is a problem of "eigenvalue type"; recalling 7| e=0 = g, we will strive 
to choose 7 = j(e,g) in a (^-dependent neighbourhood, say 

h-g\<g/2, (3.2) 

of its unperturbed value g, such that (13.11) has a nontrivial solution. That we succeed is 
the content of Theorem 13.31 Consequently, our 7 will depend analytically on e, nicely 
controlled by (7 — g\ < Cg\e\. 

The subtlety of proving Theorem l3.3l lies in solving a small denominator problem. We 
go about dealing with it using a Renormalization Group method, treating such small 
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denominators scale by scale. Here we show that the framework of [BGK99] is applicable. 
The proof, though, is self-contained. 

First, view the map X \— > Q(X) as the map that takes the pair ($, ty) to 
(fi$($,\&),f^($,\l/)) witri the components = g 2 sm§ + \d <f> f($,6 + and 

Q w (<3>, vp) = A d^ f(®, 6 + \&). Then the component form of (|3.1|) reads 

ft) it) ft)- (33) 

In each entry, stands for the matrix (dbd a f)($o, 9 + ^o)- 
From (13.31) we get for the equation 

= [(P + 7 ) 2 - A/^] 1 (A/^$i) =: J$ 1; (3.4) 

Here J is a well-defined bounded linear operator from £>* to £>*, provided that eo is 
small. Checking this is straightforward implementation of Neumann series and the fact 
that the operator (T> + 7) -2 has the diagonal Fourier kernel 

{V + 1 y 2 {p,q)=5 m {iuj-q + 1 y\ p,qeZ d . (3.5) 

Using (13.21) . one obtains the bound 

hII/^b*;*?*) — (2-6) 

Remark 3.1. The definition of J is an instance where demanding smallness of e : = \g~ 2 
is natural, indeed necessary. 

Consequently, using (13.41) . we get for $i the equation 

[(V + 7 ) 2 - <? 2 ]$! = (7 2 (cos$ - l)$i + A/^$i + A/^JSi =: (3.7) 

Recall that $o| e= o = by Lemma I2T41 Therefore -ff| e=0 = 0, and $i| e=0 = 4 (due 
to $>°(z) = 4arctanz) is a physically motivated nontrivial solution to (13. 7p . In other 
words, the differential operator (T> + g) 2 — g 2 is singular. On the other hand, when 
e 7^ is small, we know that $o remains close to zero, making the whole right-hand 
side in (13. 7p small. We then hope to find a Lyapunov exponent 7, close to g, such that 
(D + 7) 2 — g 2 — H stays singular and the equation still admits a nontrivial solution close 
to the constant function 4. 

It follows from (I3.6P that the operator H appearing in (13 .7p . which lies in £(£>*) = 
£(£>*; £>*), has the useful properties below. The proof comprises Subsection 13.11 

Lemma 3.2. Denote the kernel of H E £(i3*) by H(p,q), (p,q) E 7L d x 7L d . For 
|$5m/t| < g/3, there exists an operator H(k) E £(£>*) related to H by 

{t s H){p, q):=H{p + s,q + s) = H{u ■ s;p,q), sE Z d . 

Let < a' < a. The kernel H(n;p, q) is analytic on 

{(K,e,g,j) I |9fm«| < g/3, (e,g) ED, \j - g\ <g/2} 

and it satisfies the bound 

\H(K;p,q)\ < Cg 2 \e\e- a '\ p ~ q \ 



CONSTRUCTION OF WHISKERS FOR THE QUASIPERIODICALLY FORCED PENDULUM 17 



with C = C(cr'). As for the k- derivatives, 

\H^{K;p,q)\ <Ck\(g/3-\SmK\)- k g 2 \e\ 2 e- a ' lp ~ ql , k > 1. 



Moreover, 



|-#(0;0,0) 



1-217-^1/(7 



3.1. Proof of Lemma 13.21 

Proof. To simplify notations, we decompose 

H = H 1 + H 2 with H 2 = \f M J. 

Let $ and \1/ be arbitrary functions in the spaces £>* and respectively. 

#i acts as ordinary multiplication: Hi§(9) = Hi(9)&(9) with Hi(0) 6 C. We write 
ifi for the Fourier transform of the map 9 \— > Hi(9). Denoting a kernel element of the 
operator Hi by Hi(p, q), we have Hi(p, q) = Hi(p — q). We gather that 

t s H x = Ht (3.8) 

holds, and that the kernel of Hi satisfies 

| J ffi(p,g)|<qA|e- ff b'-«l, p,?6Z d . 

Here a > is the width of the analyticity strip around the real T d of the map 9 t— > Hi(9), 
i.e., of X . Since, by Theorem 12.11 X is analytic with respect to (e, g) G D, so is 
Hi(p,q). 

Observe that the expression defining J in (13.41) may be cast as 

J<I> = [i - (p + 7 )- 2 (A/^)] - 1 (V + 7 )- 2 (A/ v ^$) = 5A0$, 

where 5, A, and O stand for [l — (T> + 7)~ 2 (A/^)] 1 , (£>+7) -2 , and A/^, respectively. 
Assuming each index a and 6 in f a ^ stands either for <f> or ip, the reader should bear 
in mind that f a ^ refers to the multiplication operator corresponding to the Jacobian 

matrix (dbd a f)($o, + ^o)- Its Fourier kernel reads f a ,b{p, q) = fa,b(p — q), whence the 
translation invariance 

tsfa,b = fa,b- (3.9) 

Denoting A(q) = A(q, q) = (iu ■ q + 7)~ 2 , we are interested in the kernel 

J(p,q) = J2B(p,r)A(r)0(r,q), p,qeZ d , (3.10) 

of J. We shall also need the "shifted version" of A(g), 

A(n;q) := (iu ■ q + in + >y)~ 2 , kgC. (3.11) 
It is related to A(g) by the property 

t s A(q) = A(uj- s;q). (3.12) 
Further, A(k; q) is analytic on {k, | |3fm ac| < g/3} x {7 | I7 — g\ < g/2} and satisfies 

|A(«;<z)| <36<T 2 . (3.13) 
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Equation ( 13.131) also means that the operator A(k) corresponding to the kernel in 
(13. lip belongs to C(B a ) with ||A(«)|| £ / Bct n < 36g~ 2 . Interpreting f a ^ as a multiplication 

operator, \\ f a ,b\\ c(B a ) - \\fa,b\\<r shows that B,0 E C{B a ). 

As in the case of Hi, O acts as multiplication by a real-analytic function whose 
modulus is bounded by C|A|. Thus, we estimate 

\0(p,q)\ < C\X\e- alp - ql and \A(p)0(p, q)\ < C\e\ e^" 5 '. (3.14) 

Bounding the kernel of B calls for the Neumann series 

oo 

B = J2 B k, with B k := {XAf^f . (3.15) 

k=0 

Clearly ||-B fc || £ , B n < (C\e\) k and \B k (p,q)\ < (C|e|) fc such that, by Fubini's Theorem, 

oo oo 

^(p) = ^B^(p) = Y,Y t B k (p,q)4f(q). (3.16) 

k=0 qeZ d k=0 

The expression of B k contains k — 1 products of the operator A A/^ with itself, which 
appear as convolutions in terms of Fourier transforms. Explicitly, 

B k (p, q) = X k Hp)UAp - 9i) • • • HQk-i)UAQk-i ~ q)- (3-17) 

Using the bound |A(p)/^(g)| < Cg~ 2 e -0 "' 9 ' we see that, for < a' < a, 

\B k (p,q)\ < (C<r 2 |A|)V <J >- <? l J2 e- {a - a ' ){lp - qil+ - +lqk - 1 - ql) < {C\e\f e^'^ . 

Thus, choosing e appropriately small we make the geometric series arising in (13.161) 
convergent and obtain 

\B(p,q)\,\J{p,q)\<Ce~^- q \ 
with the aid of (ETT41) in flBTTOD . Finally, 

\H 2 (p,q)\<Cg 2 \e\ 2 e-^ p - q \. (3.18) 

Exploiting (I3.9p . we compute 

t s H 2 = Xtstf^J) = X f M t a J = X f M {t 8 B){t a K)0. (3.19) 

With the aid of fl3T5D and fl3TT2|) . t s B k = X k [A(u ■ s)f^} k . Thus, {t s B k )(p, q) depends 
on s only through uj ■ s. Moreover, the dependence on u ■ s is analytic in a neighbourhood 
of the real line: Consider the shifted quantity 

B k (n;p,q) := X k ^ A ( K ', p)UAp ~ Qi) ' * * A ( K ; Qk-ijUAqk-i - q), 

qi& d 

which for k = uj ■ s becomes (t s B k )(p, q). The summand above is analytic on 

D g := {e | |e| < e } x {k | |9fm«| < g/3} x { 7 | | 7 - g\ < g/2}, 

and the sum converges uniformly, as is readily observed after recalling the bound (13.131) 
on A(«; q) and looking at the estimation of \B k (p,q)\. Thus, B k (K,;p,q) is analytic. 
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But the Neumann series Yl"kLo-Bk{ K '>P,q) also converges uniformly, making the limit 
B(n;p,q) analytic on D g . Evidently, (t s B)(p,q) = B{uj ■ s;p,q). The kernel Bk(n;p,q) 
defines an operator B(k). Motivated by equation (I3.19p . we extend the definition of H 2 
and set H 2 (n) := Xf^^B^A^O. Using (13. 13j) . a straightforward computation shows 
that also H 2 (K]p,q) obeys (13.181) and is analytic on D g . Furthermore, 

(t s H 2 )(p,q) = H 2 (uj- s;p,q). 

Recalling the translation invariance (13.81) of Hi, we simply take H(k) := Hi + H 2 (k). 

The bound on the derivative H( k \n;p, q) is achieved by a Cauchy estimate. To that 
end, one observes H'(k) = H 2 (k) and uses the bound (I3.18P on D g . Similarly, because 
X is independent of 7, dH/d'j = dH 2 jd^, and we get the bound on dH(0; 0, 0)/d , y. 

The constants above are independent of g, as long as < g < go. That is to say, 
the estimates hold on Uo<g<s A? = {(^^9,l) | |9fm«;| < g/3, (e,g) G D, (7 - g\ < 
9/2}. □ 

3.2. Linearized invariant manifolds: rudiments of renormalization. We now 

proceed to stating the main theorem of this section, discussing the linearization X\. 
Our proof is based on a Renormalization Group (RG) technique we present below. 

Theorem 3.3. Under the assumptions of Theorem [21 there exist a number eo and a 
map 7 = 7(e, g) on D, analytic in t, with I7 — g\ < Cg\e\, such that equation (13. ip has 
a nontrivial solution X 1 which is 

(1) analytic in \e\ < eo and 

(2) analytic in 9 in a complex neighbourhood IA of T d , 
and satisfies the physical constraint 

*iU = 4=(*i)- 
Furthermore, it is real-valued if e and 9 are real, and 

sup (0)| < C\e\ and sup |$ x (0) - 4| < Cg\e\. 
eeu eeu 

The map 7 is independent of (^o)- If X\ an d X[ correspond to Xq and Xq of Theo- 
rem \2.1\ respectively, with (^0) = and (ty' ) — (3 G M. d , then 

X'^^X^ + P). 

Remark 3.4. We chose the normalization 4, because X°(z,9) = (4arctanz,0) is the 
unperturbed solution (separatrix) and arctanz = z + 0(z 3 ). 

Remark 3.5. The pair (7, X\) of Theorem 13.31 is unique in the sense that it is the only 
one making our construction work, which is manifested by Lemma [3.111 b elow . We do 
not prove the uniqueness of X\. However, for a given solution X\ the value of 7 is 
unique: If 7' is another one, (13.11) yields {T> + 7) 2 Xi = (V + •y'YXi because Dfl(X ) is 
independent of 7. This shows that 7' = 7, because $i(0) =4^0. 

Let us commence sketching the backbone of Theorem 13.31 by recalling equation (13. 7p . 
We expand the square on the left-hand side and obtain 

(V 2 + 2jV)$i = (H + g 2 - 7 2 )$i. (3.20) 
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For a small e 7^ 0, $1 should remain close to the unperturbed value 4 = <9 Z $°(0). Due 
to the linearity of (13.201) such a solution may be normalized as ($1) = 4. Thus, we set 

$ x (0) = 4 + £(0), (3.21) 

where we demand the function £ : T d — ► R to vanish on the average, ie., 

1(0) = 0. (3.22) 

Plugging (ET2Ij) into (jCTjl results in 

(£> 2 + 2 7 £>)£ = 7r (£ + 4), where tt := H + # 2 - 7 2 - 
After switching into Fourier representation, this reads 

i(q) = G(q) [J2 MQ,P)i{p) + Po(q)] if q e Z d \ {0}, (3.23) 
= 5>o(0,p)|(p)+po(0), ( 3 - 24 ) 



where po is a function defined through its Fourier transform by setting 

p (q) :=4vr (g,0). (3.25) 

The symbol G(q) stands for the diagonal element G(q, q) of the operator G whose Fourier 
kernel is given by 

Gfe 9 )==v.{r w-, " ( "' ,),rl l q q lT {0} ' (3 - 26) 

The matter of the fact is that, in terms of our new notations, any solution £ of 

£ = G(7r £ + A)) (3.27) 

also solves (13.231) ; only the zero mode constraint f!3.22|) has been included here. After 
finding such a £, we go on to show that it is a solution to (13.241) . as well. 

As is apparent from the definition of G, this problem involves arbitrarily small de- 
nominators io ■ q. Our strategy is to recursively decompose G into parts, each of which 
corresponds to denominators up to a given order of magnitude. We then end up solving 
"partial problems" of (I3.27P scale by scale, and show that these solutions converge to a 
true solution of (13.271) as the recursion proceeds and smaller and smaller denominators 
become dealt with. 

Leaving the all-important scaling parameter N e ]0, 1[ to be decided latei0, we shall 
need the entire functions 



if n = 0. 



Their importance lies in the fact that the sequence (% n — Xn+i)neN of functions is an 
analytic partition of unity on R \ {0}; on this set < 1 — xn S 1 pointwise, as 
iV —>■ 00. Some of the first members of the sequence appear plotted in Figure [3j The 
number 6 in the exponent is a choice of convenience; it is the one used in [BGK99]. 



Aleph, H, is the first letter in the Hebrew alphabet. 
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Figure 3. Graphs of Xn - Xn+i with n = 0,1,2,3, and N 
maxima are located roughly at K n . 



\. The 



Let us now introduce the diagonal operators G n and r n , n G N, defined by 

G n {q) ■= Xn{u ■ q)G{q) and T n := G n - G n+1 , 

respectively. Observe that Go = G and G n (0) = 0. The point here is that in T n (q) the 
functions Xn{^ • q) — Xn+ii^ • <?) act as cutoffs for the values of uj ■ q. Each T n deals with 
the denominators uj ■ q that are roughly of order ^ n and, intuitively, 

n-l 

T<n : = 2_j = G — G n (3.28) 



fc=0 



gets closer and closer to G as n tends to infinity. Instead of the full equation (13.27j) . 
consider the easier, approximate problem 

x n = T <n (rr x n + p ), (3.29) 

obtained by replacing G with T <n . It is easier since T <n discards the most dangerous 
ones of the small denominators. However, its solution should become a better and better 
approximation of the solution of (13.271) with increasing n. 

Having G = G\ + T , we decompose £ = £i + r] and assume that f] = ?7o(6i) solves 
the "large denominator problem" 

f?o = r (7r (£i + ?7o) +Pq). (3.30) 
Then, solving the original problem ( 13. 27ft for £ amounts to solving 

ei = Gi(7r (ei+»?o)+A)) (3-31) 

for Ci- 

Assuming 1 — Tovro is invertibldj, we can extract 770 out of ( 13.301) and get 

770 = (1 - r 7r )- 1 r (7ro£i + p ). (3.32) 
Therefore, (13.311) transforms into 

6 = G 1 (i-7r r rW 1 + po) 



3 Think of Tq as comprising only large denominators and ttq being proportional to e. 
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with the aid of the identities 

^(l-IVo)- 1 = (1 - Trol^-Vo and (1 - n T )-\ T = (1 - Tro^)- 1 - 1. 
Thus, defining the new objects 

7Ti := (1 - 7r ro) _1 7ro and p x := (1 - 7r r )~Vo, 

we obtain 

770 = r (7Ti^i + pi) 

and 

£i = Gi(7ri£i+Pi). (3.33) 

Indeed, equation (I3.33|) has 'precisely the same form as the original problem (13.27p . 
Now, relaxing the assumption that rj be a priori known, suppose we are able to solve 
(I3.33p . and take (13.321) as the definition of r] , instead. Then the solution of the full 
problem is recovered using the simple relation 

£ = 6 + vo = (1 - rvro)- 1 ^! + r oPo ). 

Owing to the aforementioned formal covariance between equations ( 13.271) and (13.331) . 
we may iterate the construction above. Thus, in general, solving 

£n+l = G n+ i(7r n+1 £ ri+1 + p n +i) (3.34) 

for £ n+ i with the definitions 

7T n+ i := (1 - 7r n r n ) _1 7r n , (3.35) 

p n+ i := (1 - 7r n r n )~Vn, (3.36) 

r] n := r n (7r n+1 £ n+ i + p n+ i), (3.37) 

produces £ n = £ n+i + rj n , or 

£ n = (1 - T^)- 1 (£n+i + T n p„) (3.38) 

for the solution of £ n = G n (n n ^ n + p n ). 

Equations f!3.38j) and (13.361) reveal through 

7T„£„ + p n = ?Tn [(1 - T n 7r n ) + T n p n+ i] + (l - H n T n )p n+ i 

the recursion invariance 

Mo + PO = 7Tl6 + Pi = ■ • ■ = 7Tn£n + Pn = ' " " (3.39) 

in our construction. 

Let us tidy up the notation by giving the definitions 

v n {y) = n n y + Pn and f n := 1 + T <n v n with r <0 = 0. (3.40) 

In particular, (13.391) takes the form v n (£ n ) = t>o(£o)- We also set 

E n (y) = (1 - r^vrj- 1 (y + I> n ) , (3.41) 

such that (I3.38P reads £ n = E n (£ n+ i), and (13.391) reduces to 

^n+i = v n o 5„. (3.42) 

The latter is a convenient way of writing v n+ \ = (1 — 7r n r n ) u n . Notice also that S n is 
formally invertible. 
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One easily verifies 

£ n = l + I> n+1 . (3.43) 

As a consequence, 

fn+l = fn o E n . (3.44) 
Since fo = 1, we have the cumulative formula 

/ n =S oEio...oS B _ 1 . (3.45) 

Hence, a similar expansion of (I3.42[) implies 

v n = Vo ° in- 
Inserting here the definition of f n , we get 

v n = v o (1 + T <n v n ). (3.46) 

Proposition 3.6. Let £ n = H n (£ n+ i). If £ n +i satisfies £ n+ i = G n+1 (ir n+1 £ n+1 + p n +i), 
then £ n satisfies £ n = G n (ir n £, n + p n ), and vice versa. 

Proof. Suppose £ n+ i = G n+1 v n+1 {£ n + x ). By G n = G n+1 + T n and 03.421). 

G n v n o E n = G n+1 v n+1 - 1 + 1 + T n v n+1 . 

But, with the aid of f!3.43p . this transforms into 

{G n v n -l)oS B = G n+1 v n+1 - t. 

As E n is invertible with = E n (£ n+1 ), the identity above proves the formal equivalence 
of the small denominator problems (I3.34p . or G n v n (£ n ) = £„, with differing n. □ 

Recalling (I3.45p . we immediately arrive at 

Corollary 3.7. If £ n = G n (n n £ n + p n ), then 

£o := fn{in) =£n + r< n w„(^„) 

solves the complete problem: £ = Go^oCo + Po)- 

Remark 3.8. The solution £ above comprises two terms having clear interpretations. 
The first term, £ n , solves the small denominator problem, namely £ n = G n (7r n ^ n + p n ), 
at the nth step. The second term, T <n v n (^ n ), on the other hand, consists of the sum 

n-1 

"<n(£n) := ^Vk^k+i) with ( W = ( S W '"°Vl)(U 
fc=0 

where rj k = rj k (^ k+ i) solves the large denominator problem r\ k = T k v k ((, k+ i + r] k ) in 
analogy with (|3.3U|). Indeed, T k v k (£ k+1 + r} k ) = T k v k (£ k ) = T k v k+1 (£ k+1 ) = %. 

Finally, we make a crucial observation. If we operate on (13.461) by T <n and set 

x n := f n (0) = T <n v n (0), (3.47) 
we solve the approximate problem f!3.29[) : 

x n = T <n (n x n + p ). 
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We shall demonstrate that the approximate solutions x n form a Cauchy sequence in a 
simple Banach space, and that their limit 

f := lim x n (3.48) 

n^oo 

solves the original equation (13.271) . 

Just to motivate the above discussion, think of an abstract map lZ n that takes 
G n ) to (7r n+ i, p n +i, G n+ x). The recursion scheme 

£ = G(-K £, + p ) ^ £l = Gl(7Ti^i + pi) S • • • £ n = G n (7T n ^ n + p„) • • • 

is called renormalization of the problem, and TZ n is the corresponding renormalization 
transformation. Then, in view of Proposition [321 it remains for one to demonstrate that 
this process "converges", in order to be able to solve the original equation £ = G(7r £ + 
p ). That is to say, one wishes that the renormalization flow of the triplet (tt , p , G ), 
("TTn, Pn,G n ) = (Ofc=o (^Oj POi Go) , in a sense tends to a fixed point (n*,p*,G*) of 
some limiting operator "T^oo = linifc^oo 7^" as n — * oo, and that the equation 

C = G*(tt*C + P*) (3.49) 

is well-defined and solvable. 

In our case G*p* = 0, such that the equation is linear and possesses the trivial solution 
£* = 0. Corollary 13 . 71 then throws light on why (I3.48j) should solve (13.271) : /«(£«) solves 
it, and £ n approaches zero with increasing n. Therefore, it is fair to expect that also 
lim^oo fn{0) is a solution. 

3.3. Banach spaces. Technically speaking, we need to control the renormalization flow 
(I3.35p - (I3.37I) by estimating the kernel elements of r„ and n n , for the operators 1 — 7r n T n 
and 1 — r n 7r„ had better be invertible between suitable spaces. Such Banach spaces will 
be defined in this subsection. 

We begin by analyzing the properties of the operators T n . A priori, one expects the 
most significant contribution to arise from such g's that u-q = 0(R n ), due to the cutoff 
Xn — Xn+i in the definition of these operators. Therefore, (I3.26P implies 

\T n (q)\=0(g- 1 X- n ). (3.50) 

More accurately, it is fairly easy to obtain 

it n > 1, (3 51) 
1 if n = 0, 

for I = 0, 1, . . . , 6, in a strip |ShtiAc| < C only depends on the parameter b. Pay 
attention to the fact that G(q), which was defined in f )3.26p . only depends on q through 
uj ■ q. Therefore, it is handy to introduce the analytic function 

i : C \ {0, 2i 7 } C : l{k) = (2i 7 K - re 2 ) -1 . 

In particular, l(uj ■ q) = G(q) for q ^ 0. This motivates the further definition 

\X n (w ■ q + k) - Xn+i{^ ■ q + k)]l(u • q + k) ifu-q + K^O, 
if u) ■ q + k = 0. 



r„(«;p,g) : = S. 
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The importance of the resulting operator r n (/c) is based on the possibility of viewing 
u ■ q as a complex "variable" : 

T n (q, q) = T n (0; q, q) = T n (uj ■ q; 0, 0). 

We shall often write T n (n; q) instead of the complete T n {n] q, q). 

Imposing the condition b < g/2 on b we get within |3fal/c| < K n & that 

\T n {n-q ) = N — — — \{cj ■ q + k) i(u ■ q + k)\ 

|N n (u ■ q + K)\ 

<C , r ^ 1 K~ n min(l,|N- n (wg + fi:)| 5 )^ e H n > 1, ^ ^ 

[ 1 if n = 0, 

making use of (13.511) . In particular, we have confirmed the heuristic estimate (13.501) . 

Now to the spaces promised. Ultimately the solution of (I3.27p . namely £ (and there- 
fore $i) will live in the space B% C £ 1 (Z d ;C) for a sufficiently small width a* of the 
analyticity strip — see Subsection 12.11 The following weights will come in handy: 

( H~ n \u-q\ :f > 1 

»»(?):= ' "n ^ 

II if n = 0. 

We extend these to negative indices by setting W- n (q) = w n (g) _1 . 
Definition 3.9 (Spaces h n ). For n G Z, let 



These norms induce the Banach spaces ft,„. Observe that /i is the space £ 1 (Z d ; C). 
Notice that our weights satisfy 

wWi(<?) K = w n (q) and u> n (g) > 1 (n > 1). (3.54) 
The spaces at hand thus realize the embedding hierarchy 

h n+ i C h n (n e Z) 

due to the trivial inequalities 

IHI n <HU (neZ). (3.55) 

Operator norms || ■ \\c{h n -h m ) between such spaces h n and h m will be denoted by || ■ || 
for short. We actually have, 

\\ L \\n;m = SU P^2\ L (P^)\ W m(p)w^ n (q) (TO,7lGZ). (3.56) 



Either from this or from (I3.55P by the Schwarz inequality one infers 



l,i < II ■ II < II • II , 1 (m, n £ 1*) 

ln+1-.in — II Un;m — II lln;m+l V ' ^- / 



such that the operator spaces satisfy 

£{h n ; hm+i) C C(h n ; h m ) c £{h n+1 ; h m ) (m, n G Z). 
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Moreover, the Schwarz inequality implies the useful bounds 

\\LiL 2 \\ n . m < II^II^II^H^ (l,m,n G Z). (3.57) 

From now on, n will always assume nonnegative values. Define the domain 

D n := {k G C| \k\ < K n 6}, (3.58) 

recalling that b < g/2. Then (13.521) easily validates the bounds 

lir^lL^^Cr^-" (3.59) 

for k G D n , where the (new) constant C-p is independent of k and g as long as g < g . 
This shows, in particular, that 

r n («) G C(h- n ; h n ) C C(h ; h ). 

Remark 3.10. The weights w n (q) arise as follows. The diagonal kernel of T n is strongly 
concentrated around small denominators to • q of order K n ; for large u ■ q the value of 
T n (q) is very close to zero, but not quite equal to zero. Therefore, in an expression such 
as T n £(g) = r„(g)£(g) we cannot let \i(q)\ be arbitrarily large for large values of u • q. 
This "tail" can be of the order of w n (q) = e K "I"' 9 ', say, which amounts to £ G h- n . 

It has to be emphasized that having the same power of N~ n and \u ■ q\ in w n (q) is 
crucial, which can be read off from (I3.52p . This way to ■ q "scales" as W 1 in all estimates 
in the nth step of the iteration. 

The motivation for introducing the spaces h n , on the other hand, comes from the fact 
that in the recursion (I3.35|) the domain of 7r n will shrink. So, in the norms || • || n we 
incorporate a weight that increases as n grows. It is a matter of convenience to use the 
inverse of the weight uv^g)" 1 appearing in || ■ ||_ n . 

3.4. Renormalization made rigorous: estimates and the Lyapunov exponent. 

The rest of this section is devoted to demonstrating that the renormalization flow of n n 
in (I3.35P is controlled in the norms || • \\ n ._ n such that the products ||7r n || n ._ n ||r n ||_ n . n are 
small, so as to make the recursion formulae (I3.35p - (l3.37p well-defined through Neumann 
series. Recalling (13.591) . the task roughly amounts to making sure that ||7r n || decays 
at least as rapidly as K n with increasing n. 

According to Lemma [3.21 tt = H + g 2 — 7 2 G C(B®) can be written as 

tto(p, q) =Po(w- q)S P , q + vr (p, q), 

where ttq vanishes on the diagonal, and in the first term 

Po(«) : = S Q +Pq(k), po(0) = 0, 

depends analytically on k, as long as |Chti«| < g/3; explicitly S 
and Pq(k) = H(k; 0, 0) - H{0; 0, 0). 

Similarly, we split 7r n into its diagonal and off-diagonal parts: 

7Tn(p, q) =Pn{u- q)5 p , q + 7T n (p, q), n n (q, q) = 0, 

with 

Pn(«) = Sn Pn(0) = 0- 



= J ff(0;0,0) + ( 7 2 - 7 2 
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The possibility of doing this follows from the computation 

t s 7r = t s H + g 2 - 7 2 = H(u ■ s) + g 2 -j 2 = : tt (u ■ s) 
and its recursive consequence 

t s n n+1 = (1 - 7i n (u ■ s) T n (u ■ s))" 1 ir n (u ■ s) = : 7r„ +1 (^ • s). 
Motivated by the computation above, let us inductively define the maps 

n n +iA K ) := i 1 ~ ^(k) r n (K))~ 1 7r n/3 (fi:), neD n , |$m/3| < a n , (3.60) 
starting at 

tto/3(k) := P Q («) + 7f ^(Ac), re G D , \3m (3\ < a , 
by setting b < g/3 in (13.581) . Here 

P (n;p,q) := p (re + lu ■ q)6 m , 

M«5P>?) : = e^-«>#(re;p,g)(l-5 M ), 



and, with a' coming from Lemma 13. 2\ 

4 

(n + 3) 

In particular, Eric Weisstein's World of Mathematics [Wei] tells us that 



a n+ i := ( 1 - j- ■ oN2 ) a n , a < a'. (3.61' 



a ■ J] ( 1 " 4 ) = ? > as n °°- ( 3 - 62 ) 



fc=3 



As far as notation is concerned, we may omit /3 if it equals zero: vr n (/t) = 7r„o(re), and 
so forth. By a straightforward induction argument, 

7T nP (K;p,q) := e lHp - q) 7i n (n;p, q), 

such that /3 does not enter the diagonal of ir n p. Of course, 

M^,P,q)\ = e-^< p '^\n n (K;p,q)\. (3.63) 

For clarity, set 

P„(re) := 5 n l + P n (K) with P n (K;p,q) :=p n (K + uj ■ q)5 p , q , 
so that we may express the operator 7r n ^(re) itself, without reference to its kernel, as 

7l nf3 (K) = P n (K) + TT n p{K) = 5 n + P n {n) + 7t n p{K), 5 n = 5 n t, 

for short. This decomposition satisfies 

ll^OOIU-n < N + ll^(«)|| n ._ n + \\Znp{K)\\ n ._ n . (3.64) 

It will turn out that the sum in (I3.64j) is finite if re G D n and |Qtn/3| < a n — indeed very 
small, as we are trying to prove — meaning that n n p(n) G C(h n ; h- n ). 

The crux of analyzing the renormalization flow is the following lemma, for which we 
provide an inductive proof later on in this section. The reader is advised to take the 
result as granted for now. 
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Lemma 3.11 (Modified Lyapunov exponent controls the flow). Set b = g/3 and K 

> ; C > 0, c > 0, [A > 1, and a unique Lyapun 

7 - g\ < c y \e\g (3.65) 



min(|,6 2 ). There exist constants c 7 > ; C > 0, c > 0, \i > 1, and a unique Lyapunov 
exponent 7 satisfying 



such that, for any n G N, the bounds 

ll*«isWIL-„ < ckM L _ m » "'"Z? i-w-i 




IAWIU-.<C|«|, J. im (3.67) 



5 ifn = 0, 
K 2n ifn>l, 



K\<C\e\ g r 2n (3.68) 



/io/d irue /or (e, 0) GD, kG -D n and |Sm/3| < a n . Moreover, c is bounded away from 
zero and fi —>■ 00 m t/ie /zmzi — >• 0. 

Remark 3.12. The sole purpose of introducing the complex variable k is to go about 
proving the bound (I3.67H on the diagonal part of 7r n . We use analyticity in k and restrict 
the latter to a domain of ever decreasing size. 

The possibility of including the complex parameter j3 in the analysis, on the other 
hand, facilitates proving exponential decay of n n (n;p,q) in the quantity \p — q\. This 
is sufficiently rapid for obtaining the bound (13.661) on the off-diagonal part of n n . Also 
the analyticity strip of (3 around R is taken narrower and narrower upon iteration, but 
no narrower than a certain limit («o/6). 

Corollary 3.13. The bounds of Lemma \3.11\ imply 

11 / mi ^ n\ 1 J 9 i f n = 
\K^)\L. n <C\e\g^ n , fn ^ 

The caveat to get around in the proof of Lemma 13.111 is that 5 n is reluctant to go 
to zero along the recursion. To change the state of affairs, we fine-tune the Lyapunov 
exponent 7 such that also 5 n — > as n — ► 00. As stated in the lemma, there turns 
out to exist precisely one such value of 7. This is what ultimately enables us to prove 
the convergence of our renormalization scheme, consequently validating Theorem 13.31 
discussing the linearized solution X\. For the sake of continuity, we first give the simple 
proof of Theorem 13.31 and only then prove Lemma 13.111 

3.5. Proof of Theorem E23 With oCfi cts in (13.471) . the task is to show that the limiting 
function £ — see (13.481) — is an analytic solution to (13.271) . 
Given the formal definition yp := r^y, (I3.47P implies 

X n /3 = f 71/3(0). 

Recalling (I3.44p and (13.451) . one clearly has 

fn+1,/3 = fn/3 ° ^n/3 and f nj 3 = S 0/ 3 O ^ lj3 O ■ ■ ■ O r, n _ 1( g. 
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Hence, the recursion relation 

x n +i,p = x n /3 + (/„/3(H n/3 (0)) - / n/ a(0)) 

follows. Here (13.411) extends to 

E n p(y) = (1 - T n n n p)-\y + T n p nP ). (3.69) 

Notice that the flows of p n( g and 47r nj g( • , 0) are identical. Furthermore, the initial 
conditions agree according to (I3.25p . such that 

Pnp(q) = 47r n/3 (g,0). 

Because DE, n p{y) = (1 — r n 7r ni a) _1 , the chain rule reveals 

Df nP (y) = (1 - ToTTop)' 1 ^ - Tm^)- 1 •■•(!- T^tt^)- 1 . 

Recursive implementation of Corollary 13.131 in the form 

||(i - rvr^r 1 !!^ < ||(i - r n ir n p)-%. n < 2 

implies that Df n/3 (y) G C{h n -i,h Q ) with $nv yehn \\D f nP {y)\\ n _ VQ < 2 n . By the Mean- 
Value Theorem we go on to estimate 

IK+i i/3 - ar^Ho < 2 n ||~ n/3 (0)|| n (3.70) 

with the aid of the inequality || ■ < || ■ || n . 

Lemma 3.14. For parameters as in Lemma \3.11\ and eo small, we may perceive as 
an analytic map from h n to h n C fo„_i with 



^(0)|| n <C|e| 



9 if n = 0, 
ifn>\. 



Proof. Since T n annihilates the zero mode (r n (0) = 0), T n p n p = 4(r„7f n| g)( ■ , 0), which is 
super-exponentially small in the norm || • || n by ||7r n /3( ■ ; 0)||_ n < ||7r n /3|| n _ n and Lemma 
13. Hi According to ( 13.691) . Lemma [3.141 clearly holds if we take e small enough so as to 
validate HrjJ^.JlTTn/gH^.^ < |, say, for each n. For the bounds on Tr n/3 and r„ we refer 
the reader to Corollary 13. 131 and (I3.59p . respectively. □ 

By Lemma [3.14[ f n p maps h n _\ to ho, confirming that x n p G h for each n. Coming 
back to (I3.70P and taking |Sftti/3| < a* := cto/6 (see (13.621) ). the sequence (x n ,3)neN is 
Cauchy in the Banach space ho. Moreover, xop = gives us 

oo 

ll^llo < W Xn +^P ~ x "P\\o ^ c \ e \9, 

n=0 

where the factor g is due to the last statement in Lemma 13.111 It is implied that 

Uq)\ < C\e\g e^ (qeZ d ). 
We infer that £ is real- analytic on T . 



7Tn/3( • 1 0) is shorthand for the function 9 i— > ^ e lq ' e ir n 0(q, 0). 
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Recalling that lim^oo T <n (q) = G(q) for each q G Z d , let us take the pointwise limit 
n — ► oo of (I3.27P in the Fourier representation: £(g) equals 

lim T <n (g) (5S5;(g) + po(<?)) = C(<?) hm (ttqx;^) + p (<?)) = G^C^of + A>)(?)> 

n— »oo n— too 

because 7r is a continuous operator on h . Indeed, £ solves (13. 27ft ! 
Out of curiosity, we conclude by the recursion invariance (I3.39|) that 

£ n = G n (ir n £ n + p n ) = G n (7i £ + po) (3.71) 

converges to = 0, pointwise in terms of the Fourier representation. Hence, equation 
(I3.34p really trivializes in the large- n-limit. Another way of seeing this is the pointwise 

bound \G n p n (q)\ < C\G(q)\\\ir n \\ n ._ n , which tends to zero and paraphrases G*p* = 
below (jg2BD . 

We still need to demonstrate that the solution £ of (13.231) also solves (13.241) . i.e., that 
(v5 + Po)(0) = tt (0, • )| + p (0) = J2 M0,P)i(p) + Po(0) = 0. 

From dHZZH), Gn{q) ■= XnioJ ■ q)G(q), G(0) = 0, and (E2ZD, 

|Sn(0)| < ^\n n (0,q)\Xn(w-q)\G(q)(Q + po)(q)\ 

q£Z d 

< U\\ sup\iT n (0,q)\xn(uJ ■ q) 

q£Z d 

< UWo ll^n|| n; _ n sup u> n (g)x n (o; • q) 

q&Z d \{0} 

< C U\\o ll^nllft;-* " as n^OO, 

Thus, 

(ttoC + Po) (0) = lim {ji n ^ n + p n ) (0) = lim p n (0). 

n— »oo n— >oo 

But 

lim p n (0) = 4 lim 7r n (0, 0) = 4 lim 5 n = 0, 

n^oo n— >oo n— >oo 

and we are done with the construction of (7, Xi) under the assumption (^> ) — 0. 

The case (\&o) 7^ 0. If X\ solves (13.11) . it is a matter of applying the translation 773 on 
both sides of the equation to get (V + i) 2 X[ = Dn(X^)X[, where X' Q = 773X0 + (0,/3) 
and X[ = TpX]_. In other words, the translation property in the formulation of the 
theorem holds, and the value of 7 does not change under such translations. □ 

3.6. Proof of Lemma 13.111 We begin by deriving several identities that are easy to 
refer to below. To this end, let us look at the flow (I3.60P more closely, observing that 
we may formally split 

(l - -K n p(K)T n (K)Y l = (1 - P n (K)T n (K))~ x + r nf3 (K). 

The remainder r n( g(/c) reads explicitly 

r n p{K) := (1 - 7r n/3 (/t)r n (/t)) -1 7r n/3 (/t)r n (/t) (1 - P n (/t)r n (/t)) _1 . 
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In fact, this quantity is asymptotically very small in C(h- n ; h- n ) due to the explicit 
factor n n p] given the bounds (13.66 1) f |3 .681) for some particular value of n, 

IWK)||_ n; _ n <C|£|e- c ^. (3.72) 

Continuing abstractly, (13 .60 j) becomes 

7l n+ltf3 (n) = (1 - P n (K)T n (K)Y l P„(/c) + S n (K) + S n/ g(Ac), 

where s n (/«) is the diagonal and s n p(n) the off-diagonal part of the small remainder term 
r„0(/c)7r n( g(/c), respectively. Therefore, the diagonal P n (ft) — containing the problematic 
5 n — and the off-diagonal 7T n p(K) iterate according to the rules 

(1 - P„(/t)r n (fi:))~ 1 P n {n) + s„(k), 73 
= s n p{n). 

Notice that s n («) is indeed free of /?, because each P n (/c) is. 

By construction, <5 n = 7r ni g(0;0,0) = P n (0;0) for each n, such that the diagonality 
of (1 — P„(0)r n (0)) _1 with r n (0; 0) = implies that changes in 5 n upon iteration only 
arise from the small term s n in (13.731) : 

S n+ i = S n + d n , d n : = s n (0; 0). (3.74) 

But r n/3 (0; 0, 0) = 0, again because r„(0; 0) = 0, such that 

dn = s n (0; 0) = (r n vr n ) (0; 0, 0) = (r n # n )(0; 0, 0). (3.75) 

We remind the reader of our convention of dropping one of the kernel indices of diagonal 
operators. For instance, s„(k; q) = s u (k; q, q). 

It is convenient to spell out a consequence of (13. 73ft : 

P n+l {K) = P n {K) + P n (K)T n (K)(t - P n {K)T n {K))- l P n {K) + (s„(«) - d n ). (3.76) 

Proof of Lemma VJ. 11\ Here we finally prove that the bounds f)3.66l) - fl3.68l) . such that 
(I3.73P — and indeed everything above — becomes not only formally justified. To this end, 
we proceed by induction. As iterating (I3.66|) and (13.671) is rather easy, the proof boils 
down to choosing the value of our free parameter, the Lyapunov exponent 7, so as to 
guarantee that 5 n satisfies (13.681) at each step. 

(i) Case n — 0. Consider k restricted to D with b < g/3. Lemma l3~2l and Po(«; q) = 
H(k; q, q) — P(0; 0, 0) readily imply 

||Po(k)|| 0;0 <G)|A|. 

Furthermore, increasing C and employing (13.631) with l^m /3| < a < a', 

||^o/3(k)|| ;o < Co|A|. 

The leading Taylor coefficient p' Q (0) = P'(0;0,0) and the corresponding remainder of 
the function p = H( ■ ; 0, 0) - H(0; 0, 0) satisfy 

\P' (0)\ < iCokl 2 ^ and \p ( K ) - P ' (0)k\ < ^C \e\ 2 \n\ 2 , 

taking Co large enough. 



32 



MIKKO STENLUND 



Assume that 7 lies in the open ^-centered disk of radius c\e\g: 

7£V=%* (3.77) 

Recall that 5 = eg 2 u(e,g,-f)+g 2 --? 2 , where eg 2 u(e,g,-?) = #(0;0,0). If 5 (7i) = $0(72) 
and we denote 7^ = g(l + xt), the Mean- Value Theorem yields 

1 7i -7a| < i(\xi +x 2 \ + |e|flr||9 7 it|| 00 )|'yi -7 2 |. 

By Lemma 13.21 H&ywH^ < C|e|(7 _1 /(1 — 2c 7 |e|), and |xi + a; 2 | < 2c 7 |e|. For a sufficiently 
small |e|, we gather 7! =72, such that 7 1— > 5q is one-to-one on J 7 . Moreover, the image 
of the disk J 7 contains the disk D(0, (2c 7 — c 2 \e\ — IMI^JIelp 2 ). Thus, for a sufficiently 
large value of c 7 and small value of e, there exists a closed set Jo C J 7 which 7' 1 — ^ <5 
maps analytically and bijectively onto the closed disk 

I :=D(0,C |e|s 2 ). 

We are about to prove below that a correct choice of 7 leads to 

5 n E I n :=B(0,C \e\gX 2n ) 

for each and every n 6 Z + . 

(ii) Induction step: hypotheses. Fix nGN. Suppose 



g if n = 0, 
e"^" ifn>l, 



||7f^(«)|| n; _ n <C n ieb^ 
for some constants c > and fi > 1 — to be fixed later — and 

||Pn(«)|L_ n <C n |e|#N 



<7 if n = 0, 
1 ifn>l, 



hold true for |e| < e n , \$smf3\ < a n , and k G D n . Suppose there exists a closed set 
J n C I-y and a bijective analytic map A n : J n — > /„ : 7 1— > 5„. 

Further, let the kernel elements of these operators be analytic in D n and continuous 
in the closure D n . Also the estimates 

wm\< (i-r^)^<i;:L l! 72 :?' (3.78) 

and 




n + 2; 2 ny I |e| 3 / 2 ifn>l, 



1 If I 2 if n = D 

|Pn(K)-p;(0)«|<-C n ^|K| 2 ^ ' (3.79) 

e ' it n > 1 



which facilitate dealing with the Taylor expansion of p n , are supposed to be satisfied. 
In particular, it follows from (13.641) . b < g/3, and the inductive hypotheses that 

\ Pn (K)\<B n C n \e\g^ with B n :=n e \+ 9 * n = °> (3.80) 

I b\e\ 1 + W it n > 1, 

and 

lk^(/c)|| n ._ n < A n C n \e\gtt n , (3.81) 
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where 



. q if n = 0, , . 

A n :={ y „ ' 3.82 

'l + K n + e- c ^ ifn>l. V ; 



The strategy is to iterate the above hypotheses and prove that, in the bitter end, C n 
and e n can be chosen in an n-independent fashion, uniformly in g. 
(ii a) The off-diagonal jt n +i^(K). If $ G C d , then 

But with a modification of (13.721) . 

ifn>l' (3 ' 83) 

such that 

\\7r n +lA K )\\n;-n < \Vnp{^np{^)\\ n -- n < 2C n \e\gK n B n 

both provided e meets the condition 

|e| < e n+ i := max (e n , KA^Cr)" 1 ). (3.84) 

Hence, if |Qtn/3| < a n , 

K + i,p(K; p,q)\ < 2C n \e\gH n B n ■ (1 - 5 m ) e^~^ p -^w n (p)w n (q). 
Now assume |Qttl/3| < a n+ i and, fixing p and q, take 

(3 = (3 + Z(a n - «n+l)l [• 

\p-q\ 

Obviously |Shn/3| < a n . What we get this way is 



7T 



5 («;p,<z)| < 2C n \e\gX n B n ■ (1 - 5 p>q ) e^- a ^^w n {p)w n (q) 



for each pair (p, q) G Z d x Z d . Thus, from the expression (13.561) for the norm, 

ll^n+l,/3( K )lln+l;~(n+l) 

<2C n |e|^£ n sup y e -^ + 3)'^\ P - g \ Wn(p)w n (q) 

le&p&^iq} w n+1 (p)w n+1 {q) 

<2C n \e\gK n B n ]T e^^'^w^AvY^ ■ (3-85) 
pez d \{o} 

After (13.541) . the second inequality follows from shifting p to p + q. We control the above 
bound by treating the cases \uj • p\ < M;( n+1 )/ 2 and \u • p\ > ft( n+1 *>/ 2 separately. In fact, 
if \wp\< K( n+1 )/ 2 , then \p\ > ^-(^+1)/^ follows from (031) . and 

e _ 4 (n+3)- 2 a„|p| < e -2 n - 2 a n \p\ . e -2(n+l)- 2 a„N-(™+ 1 )/ 2 ^ ^ j ^ -(1-N) < ^ 

|w • p\ > K {n+1)/2 =► Wn+1 (p)- (1 - H) < e -d-«)«- ( " +1)/2 . 
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Since a n > ao/6 by f 1 3 . 6 2 j) and, for a > 1 and m > 0, m 2 a m > ^-(\na) 2 , we have 

-2(n+l)- 2 anN _(n+1)/a " < e -^e 2 a ln(K- 1 /^)K-(™+ 1 )/ 4 " 

The remaining <i-dimensional geometric series satisfies 

J2 e -2(«+D- 2 ao|p| < C{d) f^lT) \ 

P ez d \{o} V / 

Hence, we infer that if |Shti/3| < a n+ \ and k E D n , then 

\\Zn + lM\\n + l;-(n + l) < C n+1 \e\g , 

where we finally pin down the values of the previously free parameters 

c:= imin(^e 2 a ln(K- 1/4iy ),l-N) > and := K~ 1/max(4y ' 2) > 1, 
and take ^ 

CWi > 2C(d) ^ 1 e- c ^ +1 B n ( {n + 1)2 ) C n . (3.86) 

(ii b.l) The non-constant part P n+1 (/c) of the diagonal. If k E -D n +i and 
\u-q\< K n (l - N)6, then k + cu ■ q E D n . So, by (ET53jl . 

\\Pn+i(K)\\n+i:-( n +i) = sup |P„+i (k; g) | w„ +i (g)~ 2 

/ / \Pn+l(K + U ■ q)\ - , 2/N 

< max< sup —2 , sup \P n+ i(n;q)\w n {q) 1 

I |ui-g|<N n (l-*t)6 I»n+lWJ |wg|>N"(l-N)6 

< m ax{ sup IP ~"+ l(K + ;- g)l , ^-) 2 ||P„ +1 ( K )||„._„ 

But we know that the relations ||-Pn+i(«)|| n ._„ < ||-PrH-i( K )ll n — « + \^\\ and |<5 n+ i| = 
|P„+i(0;0)| < ||P„+i(0)|| n ._ n hold. Moreover^ (13TSTT) and fl3~6T?]i yield 

ll^+i(«)IU ; _„ < ||7r n+ i^(K)||„._ n < 2||7r n ^(«)|| n ._ n < 2A n C n \e\gX n , 

assuming (13.841) and k E D n D D n+ i hold. Observe that, for positive x and p, 
x -\ e -x p /(ep) < i Consequently, if we demand that 

N<min(§,6 2 ), (3.87) 

say, and 

C n +i > A n C n , (3.88) 

it remains to be proven that 

\p n+1 (K + uj ■ q)\ +1 

sup — < C n+1 \e\gi< + . (3.89) 

K&D n+l ui n+1 (qy 

|w-g|<N n (l-N)6 

Notice that the rather arbitrary ( 13. 87ft imposes an interrelation between K and g, which 
is needed in the limit b < g/3 — > 0; since we cannot take b large, we have to take 
K = o(b) in order to guarantee e~ 2m < K/4 above. 

In order to verify (I3.89p . we use the recursion formula 

Pn+l - Pn = Pn In dnPn + S n (-]0)- S„(0; 0) (3.90) 
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subject to 

a n := (1 -pnlnY 1 and 7„(k) := Y n {n\ 0), 
which is an advocate of (I3.76p . The bound (13.521) yields 

l7n(«)| < C^- 1 ^""^-^! 5 (k e D n ). (3.91) 

By virtue of |s n (/c;0)| < ||r„(/c)7r n (/c)|| n ._ n , (I3.83P gives 



|s„(k;0)| < 4^Cr|e| 2 „ lf W °' (3.92) 

1 nX ' J{ ~ n 1 1 \A n gX n e- c » ifn>l, V ' 



in £>„. 

(ii b.2) The Taylor expansion of p n+ i(n) = P„ + i(k;0). Let us abbreviate 

<T n («) =Pn(/«) -#,(())«, 

for each natural number n. The objective is to show that the estimates 

i.e., the iterate of (I3.78p . and 

sup |K +1 - <r n )(«)l < L7 n+1 | e | 7 / 4 ^K" +1 (3.94) 

K,eD„ 

hold. Indeed, with the aid of such bounds together with (13. 79ft . (I3.89f) follows from 

sup {x + I k\ ) k e~ ax = (-) (a > 0) 

for = 1, 2 and e suitably small. Moreover, the Cauchy estimate 

\\-2n 

kn+i(«OI < |<r n («)| + b~ 2 \n\ 2 sup |(cr n+ i - <r n )(C)\ ( K £ D n+i), 

implies that also (13.791) gets successfully iterated. 
The bound in (13.911) implies 

7n(0) = 7n(0) = 0, 

such that p' n+ i(0) = p' n (0) + s' n (0; 0) according to (13.901) . and hence 

1 J s n ((;0) 

Thus, resorting to (13.921) 



_ , , , c n+1 \e\ 3/2 g 



IK+l(0) ~K(0)l < ^ & SU P |Sn(K,0)| < 



2(n + 2)(n + 3)' 



if the constant C n+ \ satisfies 

C n+1 > 8(n + 2)(n + 3)c7 r r 1 A„ J B n |e| 1 / 2 C'2. (3.95) 
The bound (I3.93P now follows, assuming also C n < C n+ \. 
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We still need to demonstrate ( 13.941) . This will be provided by ( 13.901) . since then 

K i 

cr n+ i(«) - cr n (/c) = (pnl n a n p n ){K,) + s„(k; 0) - }^ s^(0; 0) — , 

1=0,1 

such that (I3T80|) . (I3T9TD and (l3T92|) yield (gH]) if 

C n+1 > AC T K-\B 2 n b b + 6A,£„)N 1/4 C£ (3.96) 

Intuition behind (ii b.l— 2). Due to the super-exponential decay of s„(k;0) in 
(I3.92p and the strong induction hypothesis \5 n \ < Co\e\gti. 2n on the constant part of p n , 
the flow of the remainder p n = p n — 5 n reads roughly 

p n+1 PS (1 -pnln)~ l Pm (3.97) 

by ( 13.731) . Hence, the a priori bound |(1— p n 7n) _1 | < 1 H-C | e | yields a sequence diverging 
in n, with very little hope of proving bounds such as (I3.67P — see (13.891) . However, the 
support of 7fc is highly concentrated on the annulus < \k\ < Iterating for 

n > k steps, with k on the latter interval, 

p n+ l(K) ~ (1 - Pl{K) lk {K))- l p X {K) = (l + 0(e))pi(«). 

That is, p n remains close to pi, which enables proving (13.671) through ( 13.891) . 

In fact, our argument is different still: since x„(Kk) = Xn-i(^) and G(N 0) ~ 
N _1 G(k;0), we have 7 n (K«) ps ^ _1 7„-i(k) for n > 2. Inserting this into (I3.97j) . we 
notice that the approximate scaling invariance 

p n+1 (ttn) ~ Np n (/«) 

is consistent with the flow. This is what the bounds ( I3.78l) -( l3.79p reflect. 

(ii c) The constant part 5 n+ \ of the diagonal. Recall that 7 may be viewed as a 
function of 5 n by the induction hypotheses; the identity 5 n = A n (7) is bijective on J n . 
The flow produces a near-identity analytic function S n+ i = 5 n + d n (5 n ) of 5 n on the disk 
I n , such that, for e small enough, 

5 n+1 {I n ) D I n+1 . (3.98) 

The analyticity of the map 5 n 1— > c? n can be read off (I3.75P and the expression of r n . As 
far as estimates are concerned, 



|4| < ||rn(0)# n (0)|L .^<C(%Cr\e\*g 



if n = 0, 

^n e -2c^ if n > I 



in the complex neighbourhood 2I n of I n of radius \\I n \i where \I n \ is the diameter of 
the disk I n . Consequently, a Cauchy estimate yields the bound 

sup \dd n /d8 n \ < SUP ^^ 2/ ; Mnl < I (3.99) 
on the Lipschitz constant of d n on /„, provided |e| < e n+ i with 



^>2C^CClC r \l^, (3.100) 
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In this case also 

Kl < Wn\ (3-101) 

holds, which validates (I3.98p . considering how the boundary of I n is transformed under 

<Wi- 

Notice that (13T99D implies 

\5 n+1 (x) - 8 n+1 (y)\ > ±\x - y\ (x,y e I n ), (3.102) 

meaning that 5 n i— > is one-to-one. By continuity and (I3.98p . there exists a closed 
set J n+ i C /„ that is bijectively and analytically mapped onto I n+ \: J n+ i '■— 5~l 1 (I n+ i). 
We can backtrack with the aid of the map A n , obtaining a closed subset J n+1 C J 7 
(see (13.771) ) that is bijectively and analytically mapped onto I n+ i by the map A rt+1 : = 
o n +i ° A n : 

Jn+l '■= A~ +1 (/ n+ l). 

It follows immediately that 

(iii) Large values of n and the limit g — ► 0. Suppose C n is independent of 
which is the case for C . The recursive conditions (I3.86p . (13.881) . (I3.95p . and (I3.96P can 
be summarized in bounds of the form 

Cn+i > K n {g) C n and C n+1 > L n (g)\e\ 1/4 C 2 n . 

Choosing b := g/3 (due to fl3T95D ) and K := min (|,6 2 ) (due to fl3~96l) : see also fl3~87]) ). 
which is allowed, we may bound K n (g) and L n (g) uniformly in g: 
sup 0<9<90 K n (g) < K n and sup 0<9<go L n (g) < L n . This follows from the fact that 
^-lg-e^t ^ o as K ^ 0. Moreover, L n < L for each n, such that we may choose 

C n+ i : = max (K n , L|e| 1/4 C n ) C„. 

The numbers K n > 1 converge to unity so fast that the number K := Yl^=o Kn > 1 
is finite. Now choose e so small that L^^^KCq < 1. In particular, C\ = K Cq, and 
inductively C n — K ■ ■ ■ K n _iC < KC . We conclude that the sequences (C n ) and (e n ) 
(see (I3.84p and (I3.100P ) converge to positive numbers. 

(iv) Fine-tuning the Lyapunov exponent 7. The maps 5 n are relatively expan- 
sive; (I3.102p holds, while the target I n contracts by a factor of K 2 < | at each step. 
Thus, demanding A n (J n ) = I n at each step for the map A n = 5 n o ■ ■ ■ o 5 amounts to 

\x-y\< 2 n \A n (x) - A n (y)\ < Cg(2K 2 ) n (x,y e J n ), 

or limn^ool J n \ = 0. Because the J n 's form an ever decreasing chain of closed disks, their 
intersection consists of precisely one point: 

00 

{7} == n Jn c ^- 

n=0 

The value of 7 is an analytic function of e, because the sequence A~ 1 (0) converges 
uniformly to 7 with respect to e. For real values of e, A n sends real numbers to real 
numbers, making 7 real. □ 
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4. Proof of Theorem CD 

Let us summarize what we have learned thus far. The solution X°(z) to the equations 
of motion in the uncoupled case was found. In the coupled case we resolved KAM-type 
small denominator issues, which contributed the t — ► — oo (z = 0) asymptotic Xq(8) of 
the general solution X(z, 9), as well as the linearization X\{6) = d z X(0, 6). 

We can now solve (11.111) . and thus find the unstable manifold also "far away" 
from the torus T\, by a Contraction Mapping argument. 

To begin with, we single out the uncoupled part X° of the complete solution X; 

X = X° + X with X| e=0 = 0. 

As C 2 X° = ( 7 2 sin $°, 0), ffTTTD now becomes C 2 X = -( 7 2 sin $°, 0) + fi(X° + X). In 
other words, the map X has to satisfy 

KX = W(X), (4.1) 

where we define the linear operator 



^ := (o C 2 ) With ^ : =^ 2 -7 2 cos$° 



(4.2) 



h k (6) := d * h ^ 6 \ h< k (z,9):=Y,z k h(0) and S k h := h - h< k ^. 



and the nonlinear operator W through the expression 

W{X) := (- 7 2 sin $° - 7 2 (cos $°)$, 0) + fi(X° + X). (4.3) 

Throughout the rest of the work, we shall refer to different parts of the Taylor ex- 
pansion of a suitable function h(z, 6) around z = using the notation 

k\ 

Observe that X = X and X\ = (4, 0) + X\ exist. Setting 

X(z, 6) = X< x {z, 6) - (4, 0)z + Z(z, 9), (4.4) 
we may transform equation (14.11) into the equation 

K.Z = W(Z) (4.5) 
for Z = 52X, where we define W through 

^(cos^ )^^!" 



W(Z) := 5 2 



W(X) 







(4.6) 



taking now (14.41) as the definition of X. 

Let us consider the complex Banach space A of (bounded) analytic functions Z on 
the compact set 

n r := { (z, 9) (z, 9) e [-1 - r, 1 + r] x T d , 3m (z, 9) e [-r, r] d+1 } , 
r > 0, equipped with the supremum norm, and its closed subspace 

Ai := {Z G A I Z<x = 0}. (4.7) 
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For future use, let us also define the closed origin-centered balls 

B(R) := {Z e A\ \\Z\\^ < R} and B l {R) := B(R)nAi. 

Any element of A extends analytically to I1 T / for some r' > r, allowing uniform estimates 
on its derivatives on U T . 

Remark 4.1. Whereas equation fl4.ll) is plagued by small denominators, equation (14.51) 
is not. This is so due to the decomposition (14.41) which separates the previously solved 
"KAM-asymptotics" X<± from X and enables reducing (14. II) to (14. 5 p on the space Ai, 
which one could well call the small-denominator-free subspace of A. 

4.1. Existence and uniqueness of Z. Postponing the proofs until the end of this 
section, we make two observations, important in demonstrating that (14. 5p is solvable. 

Lemma 4.2. With sufficiently small R, r, and e (depending on the analyticity region 
of f), the operator W : A — > A\ maps the ball B(R) in A into a ball Bi(R') in A\ 
with R' = Cg 2 (R 2 + |e|), and W\a 1 is Lipschitz continuous on Bi(R) with a Lipschitz 
constant proportional to g 2 (R + |e|). If the restriction of Z G A to a real neighbourhood 
of [—1, 1] x T d has the real range R x M. d and e is real, then the same is true ofW(Z). 

-> A± has a bounded inverse 
It preserves analyticity in e. 
1, 1] x T d has the real range 

We have developed enough machinery to extract a solution from (14.51) : 

Theorem 4.4. For sufficiently small R, €q < R/2, and r (depending on the analyticity 
regions of f and X<\), equation H4.5[) has a unique solution Z G Bi(R). It is continuous 
on D, analytic in e, and bounded uniformly by C\e\. The restriction ^|[_i i i] X T d takes 
values mix R^, provided e is real. 

Proof. We know by Lemmata 14.21 and 14.31 that K,~ 1 W maps Bi(R) into itself. We may 
furthermore choose eo and R such that the operator /C -1 !^ becomes contractive on 
Bi(R). The Banach Fixed Point Theorem implies that K~ 1 W has a unique fixed point 
Z in the ball B^R). 

The theorem also implies that Z is analytic in e. Namely, Lemma 14.31 says that 
/C _1 preserves such a property. Furthermore, the e-dependence of W comes solely 
from 7, Xq, X±, and Q, making it analytic. Hence, the uniformly convergent sequence 
((/C _1 H^) fc (0)) fcgN reveals the analyticity of the limit Z. The latter is also R x Revalued 
on — 1, 1] x T d if e is real. Finally, 

II^IL < ||(x:-V)(z) - (xr^xoiu + W^wmiL < W\L + c\e\ 

yields H^H^ < C|e|/(1 — L). Here (/C~ 1 H / )(0) was bounded using R' of Lemma |4T^1 at 
R = 0. □ 



Lemma 4.3. If < r < 1, the linear operator fC : A\ - 
K,- 1 G C{A X ) obeying W^W^^ < C-/- 2 t-\1 - t 2 )~ 2 . 

If the restriction of Z G A to a real neighbourhood of [— 
R x R d , the same is true of fC~ l Z. 
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4.2. Putting it all together. To reach the statement of Theorem [T] about X u , we 
glue together the pieces provided by Theorems 12.11 13.31 and 14.41 

Assuming (ty Q ) = 0, we have constructed analytic maps 7 and 

X(z, 9) = X (6) + zX x {B) + S 2 X(z, 9) with 5 2 X = Z + 5 2 X° 

that solve (11.111) in a complex neighbourhood of [—1,1] X T d and satisfy the physical 
constraint $i| e=0 = 4. Recall now (11.161) . Since (11.181) is not automatically satisfied, we 
are required to pinpoint specific values of a and j3 so as to fulfill X a ^{l, 0) = (n, 0). To 
this end, we utilize the Implicit Function Theorem. 

Consider the implicit equation X(e, g\ a, (3) : = X(a, (3) + (0, f3) — (n, 0) = 0. Both 3C 
and gn^ps are continuous, and we get from X = ($°, 0) + (9(e) that 

m« 1.0) = and det (^|^) = _J_ + 0(£) 

for (e, g) G D and for whichever values of a and (3 the map X is well-defined. Hence, 
if we choose eo small enough, there exist unique continuous functions a and (3 on D, 
analytic with respect to e, such that a(0,g) = 1, (3(0, g) = 0, and 

X(t,g;a(e,g),(3(e,g)) = 0. 

Moreover, a(e,g) G M and f3(e,g) G for e real, as X is then real-valued. A good 
reference here is [Chi96]. 

4.3. Proofs of Lemmata 14.21 and 14.31 We conclude the section by presenting the 
proofs of Lemmata 14.21 and 14.31 used in the proof of Theorem 14.41 



Proof of Lerama \4-S\ Given Z G A with H^H^ < R, we study W(Z) — defined in (14. 6p . 
and clearly an element of Ai- Notice that in the relation (I4.4p . expressing X in terms 
of Z, the maps Xq and X\ were previously determined and are independent of Z. 
Furthermore, taking advantage of (14.41) and Theorems 12.11 and 13.31 we deduce 

WXW^KCdel+R). (4.8) 

With the aid of (11.101) . cast equation (14. 3 j) as 

W(X) := (0 2 sin($° + $) -7 2 sin$°-7 2 cos($ )i,0) + Afi(X + X). 

Recall that / is analytic on the strip |Qttl0|, ISmV 1 ! < V- Also, 3m$°(2:) = 0(t) on 
I1 T , when r< 1. Hence, owing to (14.81) . our function Q(X° + X) is well-defined for A 
and R sufficiently small and the strip n r about [—1,1] x T d narrow enough. 
Since sin($° + $) = sin $° + cos($°)$ + C($ 2 ), in a neighbourhood of n r 

WWiX)^ < \g 2 - 7 2 | HsincD + cos^ )^^ + Cg 2 ^ + |A| \\h(X° + X)^. 

The factor g 2 — 7 2 is the reason we chose to subtract 7 2 cos($°)$ from both sides in 
equation (14. ip . Namely, \g 2 — 7 2 | = \2g + 7 — g\ \g — 7I < Cg 2 \e\. Terms proportional to 
$ are dominated by (14.81) . Thus, for e and R small (independently of g and each other), 

||W(^)|L<C^ 2 + |e|). 
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p~4|l 

In order to obtain the Lipschitz continuity of it suffices to show that Z 

X i— > W(X) is Lipschitz, as neither (W(X))<i nor X<\ depend on Z =: 82X. To 
that end, we use the Mean Value Theorem, see [Cha85], and conclude that for some 
Z =: 82X on the line segment between two points Z' =: 82X' and Z" ='■ 82X" 

\\W(X>) - wmiL < \\DW(X)\\ \\Z' - z"\\^. 

The derivative is bounded by Cg 2 (R + |e|) given (14.81) . in particular when WZ]]^ < R. 

From its explicit expression, one immediately recognizes that W preserves the class 
of functions whose restriction to [—1, 1] x T d has the real range R x R d , if e is real. □ 

Proof of Lemma \4-3\ C maps A\ into itself, and K in H4.2[) inherits this feature. 
Let us start with the "pendulum part" of /C, and solve 

Lf = g 

resorting to the method of characteristics; we write (z, 6) = ((e 7 *, $ + ujt) in order to 
obtain an ordinary differential equation (ODE). Recalling the identity (11.71) . we see that 

(<9 2 - 7 2 cos $°(Ce 7 ')) f(Ce"\ 1? + ut) = g{^\ # + ut), (4.9) 

and our task reduces to studying L t := df — 7 2 cos $°(Ce 7 *). Since a translation in t and 
d eliminates (, we can just as well set ( = 1. 

We proceed in the Fourier language. The function / solves equation (14. 9p if and only 
if for all q G Z d the functions u(t) := e llJ ' ajt /(e 7 *, q) and v{t) := e l9 ' w *^(e 7< , q) satisfy 

L t u = v. 

Noticing that cos$°(e 7 ') = 2tanh 2 7t — 1, we see that L t has got the zero mode 

Ui(t) := (cosh7t) _1 , 

i.e., L t u\ = 0. Since L t u = is a linear second order ODE, there exists precisely one 
other zero mode U2 of L t that is linearly independent of U\. Because u\(t) ^ for any 
t G R, U2 may be found by a standard procedure: 

f dt t sinh7t 

omitting any additive constant emerging from the integral. Let us express the linear 
homogeneous equation L t u = as the first order system U = AU with U := (u,u) T 
and A(t) := ( 7 2 cos $o( e7 ^ J). Then w := ^) i s a fundamental matrix solution of the 
system (i.e., w = Aw) with detw = 1 and thus 

w -l = (** ~ U A an d W (t)w-\ S )=(* Mt)Ms) ~ 

In terms of a first order system, the complete equation L t u = v reads U = AU + V, 
V := (0,f) T . Varying constants, 

'to 



U(t) = w(t) W - L (t )U(t ) + / w-\s)V(s) ds 
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Next, we take to — > — oo. In that limit u(to) = 0(e 27 * ), such that 
u(t) = / [u2(t)ui(s) - ui(t)u 2 (s)} v(s) ds. 



Equivalent ly, 

f(e^,q)= [ K9(8ie*)§(e*e'',q)e i * tM ds (4.10) 

J — oo 

in terms of the kernel 

K*(s; z) := W* 2 (z)>V*i(;ze 7B ) - W*i(z)VM*e 7 '), 
defined (by analytic continuation) on {(s, z) G R x C | z ^ {±z, ±ie~ 7S }}, where 

W»i := IP and >V$ 2 := 7 _1 P ln+-7 _1 Q, 

and 

P( z ) := ( z 2 + l)~ l z and Q(z) := z'^z 2 - 1). (4.11) 

This is so, because VlV/(e 7i ) = Uj(t). 

In a complex strip |Qtnz| < r < 1, the inequality |z 2 + 1| > 1 — r 2 yields 

< C(l -r 2 )-^-^ 71 " 1 , s<0. (4.12) 

Since /(0, g) = g(0, q) = 0, we find that (14. 10p remains true if replaces e 7 . Inserting 
all this into the Fourier series of f(z, 6) leads to 

f° ~ 

f{z,6)= K^(s;z)g(ze^ s ,9 + us)ds, (z, 9) G [—1, 1] x T d . (4.13) 

J — oo 

Here Fubini's Theorem was used, taking advantage of the bound (14.121) . Indeed, we 
may express g{z,8) = z 2 h(z,8), where h is analytic in the same region as g. Since 
\h(z,q)\ < sup n , \h(z, 6 l )|e~ T '' 9 ' < Ce~ T '^ for some t' > r, we have on IL- that 



f \K^ S ;z)g( Z e^^q)e i ^ e+ ^\ds<C(l-T 2 )- 2 7 ~ 2 \z\ 2 £ e -(-'-)H < 



OO. 



Following the line of reasoning above, solving the "rotator part" 

C 2 f = g 

amounts to integrating u — v and results in an expression like (14.131) with the kernel 

K 9 (s;z) := W* 2 (z)W* 1 (^e 7S ) - W m {z)W m {ze> 8 ) = -s, 

introducing 

Wm '■= 1 and >V* 2 := I' 1 In . 
For each index n G N U {oo} define now 

I n (z,6):= [ K{s;z) Z{ze~< s ,8 + us) ds with K := ( K * ~Y 

J-n V 
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where (z, 9) G HV and Z G A\ are arbitrary. Also denote 

KZ := I n 



1 oo • 



Since the integrand here is an analytic function of (z, 6) on the compact region Il r and 
continuous in s G [— n, 0], it follows from an exercise in function theory that I n with 
n < oo is analytic on HV; see p. 123 of [Ahl66]. As an element of Ai, Z(z,6) has the 
representation z 2 Z(z,6), where Z is analytic on il r . Accordingly, (14.121) implies 

e" 7S \Z(ze^ s ,9 + ujs)\ds < C\z\ 2 ^-^, (4.14) 

7 e^ n 

showing that I n — > KZ uniformly on IT T as n — > oo. Hence, also KZ is analytic on the 
latter region. Moreover, I n (z,8) = 0(z 2 ) as z — > 0, which by virtue of (14.141) yields 

K:Ai-> Ai. 

We showed above that if Z G Ai and KZ = Z' (thus Z' G ^4i), then Z = KZ' holds 
on [—1, 1] x T d C II T . But each side of the latter equation is analytic on I1 T and hence 
agree there, meaning that K is the left inverse of /C: /C/C = l^. A direct computation 
shows that it is also the right inverse. In other words, 

K, = K,^ 1 on A\. 

K(s;z) G K, provided zGl Thus, should the restriction Z\[_ ltl ] x jd be real-valued, 
so is (/C- 1 Z)| [ _ 1>1]xTd . 

The integrals 7„ also depend analytically on 7. Thus, according to Theorem 13.31 
they are analytic functions on the domain |e| < e . Since I7 — g\ < Cg\e\, the trivia 
7 > \g > and (I4.14p guarantee that the convergence I n — > K~ X Z takes place uniformly 
on compact subsets of D defined in (11.151) (g bounded away from zero). 

It remains to be checked that /C _1 is bounded. For Z G Ai, Z(z, 6) = J2T=2 F Zk(&) z k 
converges in the disk D(0, r) := {z G C | \z\ < r}. Using the Cauchy inequalities 
\Zu(9)\ < k\ r~ k IIZII we deduce the bound 



|^,^)|<2(|z|/r) 2 ||Z|| oo if zGD(0,r/2) 

In Il r , \z\ < R for a certain R = 1 + 0(t), such that 2;e 7S G D(0, r/2) whenever s < 
S := — 7" 1 ln(2i?/r)). The bound (14. 12ft for Kq, applies to Kq, as well. Summarizing, 

lUC-ZII < C||Z11 - jre-~d. + f S f' R \*s) < 



'""ili-rTUs J-„t 2 (i + t) 2 y- 7 Mi-T 2 ) 2 

which finishes the proof. □ 
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5. Analytic Continuation of the Solution 

Here we present the proof of Theorem [2J In the notation of Section HI the existing 
map Z = 82X solves (14.51) and, by virtue of Ws analyticity, admits the representation 



SoX = K~ l 8 7 



7 2 cos$° 




)oo 
*<i + X> (fc) fef fc 
k=0 

oo 

+ AT 1 [>> (X< x + 8 2 X) ® k - «,« (X<0 ( 



+ 



(5.1) 



k=l 

on the set II T , taking e small enough, and denoting 



w {k) l_ D k w r^ 

k\ 



and a repeated argument of such a symmetric fc-linear operator by 

(*)•*:=(* .,*), 

k times 

for the sake of brevity. Observe that we have omitted a 82 in front of the square brackets 
on the second line of (15. ip as redundant. 

Equation (15.11) may be viewed as a recursion relation for 82X. It is crucial that 

w®\ = 0(eg 2 ), (5.3) 

when (e, g) G D; see (|1.15p . Namely, any given order 82X 1 in the convergent expansion 

00 

8 2 X = Y,^S 2 X i 

i=i 

is then completely determined by X<\ and the lower orders 82X 1 (1 < / < £ — 1) 
through the right-hand side of (15. ip . Moreover, since X<\ = 0(e), only finitely many 
terms in the sum over the index k are involved. Together these facts imply that only 
finitely many recursive steps using (15.11) are needed to completely describe any given 
order 8%X in terms of X<i alone and that, at each such step, only finitely many terms 
from the fc-sum contribute. 

It is important to understand that X<i is a predetermined function. As we shall 
see, the recursion procedure will then provide the analytic continuation of each X u ' e = 
X e Kl + 8 2 X l {£ > 1) to the large region V T ^ x {\Qm6\ < a} of Theorem [3 



5.1. Tree expansion. We next give a pictorial representation of the above recursion. 
It involves tree diagrams similar to those of Gallavotti, et al. (see, e.g., [Gal94b,CG94]), 
with one difference: there will be no resummations nor cancellations, as the expansion 
in (15.11) contains no resonances and is instead well converging. This so-called tree 
expansion is needed for bookkeeping and pedagogical purposes; we simply choose to 
draw a tree instead of spelling out a formula. 
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Let us first define the auxiliary functions 

V°)+ [( 72C0 S *°°) +W«]X<! if k=l, 



y*>(Xa) 

and make the identifications 



-0 :=K- 1 8 2 h {k) and 



Furthermore, let 
and, for k > 1, 



if A; = 2, 3, ... , 



(5.4) 



fc=0 




fc lines 



In the diagram representing the fc-linear w^ k \ the "free" lines to the right of the node 
stand for the arguments. We say that these lines enter the internal node, whereas the 
single line to the left of the node leaves it. For instance, 



O = KrV 3 )(x< l5 8 2 X, K-%h®). 




Notice that, as w' ' is symmetric, permuting the lines entering a node does not change 
the resulting function. We emphasize that all of the functions introduced above are 
analytic on Il r and |e| < eo- 

In terms of such tree diagrams, or simply trees, equation (15. ip reads 




using multilinearity to split the sums X<\ + 8%X into pieces. Above, the sum after 
the first tree consists of all trees having one internal node and an arbitrary number 
of end nodes, at least one of which, however, is a white circle. This rule encodes the 
fact that on the second line of (15. ip the summation starts from k = 1 and that the 
contributions with only X<i in the argument (i.e., trees with only black dots as end 
nodes) are cancelled. 

Using (15.11) recursively now amounts to replacing each of the lines with a white-circled 
end node by the complete expansion of such a tree above. This is to be understood 
additively, so that replacing one end node, together with the line leaving it, by a sum 
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of two trees results in a sum of two new trees. For example, such a replacement in the 
third tree on the right-hand side of (15.51) by the first two trees gives the sum 

-<!*- <r° 

Before proceeding, we introduce a little bit of terminology. The leftmost line in a tree 
is called the root line, whereas the node it leaves (i.e., the uniquely defined leftmost 
node) is called the root. A line leaving a node v and entering a node v' can always be 
interpreted as the root line of a subtree, the maximal tree consisting of lines and nodes 
in the original tree with v as its root. We call v a (not necessarily unique) successor of 
v', whereas v' is the unique predecessor of v. 

The recursion (15.51) can be repeated on a given tree if it has at least one white circle 
left. Otherwise, the tree in question must satisfy 

(Rl') The tree has only filled circles (@) and black dots (•) as its end nodes, 
together with 

(R2') Any internal node has an entering (line that is the root line of a) subtree con- 
taining at least one filled circle as an end node. 
After all, the recursion can only stop by replacing an existing white circle with a filled 
one. Continuing ad infinitum yields the expansion 

O = E ( Trees satisf y in g ( R1 and ( R2 ')) = E' t ' ( 5 - 6 ) 

trees T 

where the prime restricts the summation to trees T satisfying (Rl') and (R2'). We point 
out that each admissible tree appears precisely once in this sum, considering different 
two trees that can be superposed by a (nontrivial) permutation of subtrees that enter 
the same node. 

The earlier discussion concerning the description of 82X in terms of a finite sum 
involving only X<\ translates to the language of trees in a straightforward fashion. 
First, the second part of (15.31) and X<i = 0(e) amount pictorially to 

. =C(e) and . = 0(e). 

Second, = 0(g 2 ) and the first part of (15.31) yield 

© = (e k ) (k = l,2,...) 

and 

= 0(1) (k = 2,3,...). 
Expanding the filled end nodes 

— = E — © > ( 5 - 7 ) 

fc=i 

according to (15.41) . on the right-hand side of (15. 6ft . we get a new version of the latter by 
replacing the rules (Rl') and (R2'), respectively, with 
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(Rl) The tree has only numbered circles ((aT) with arbitrary values of k) and black 
dots (•) as its end nodes, 

and 

(R2) Any internal node has an entering (line that is the root line of a) subtree con- 
taining at least one numbered circle as an end node. 

Let us define the degree of a tree as the positive integer 

oo 

degT := #( . ) + #( .) + *#( ©) (5-8) 

k=i 

for any tree T satisfying (Rl) and (R2). By #(GQ we mean the number of occurrences 
of the graph G in the tree T. That is, the degree of a tree is the number of its end 
nodes with suitable weights plus the number of nodes with precisely one entering line. 
Since a tree has finitely many nodes, its degree is well-defined. Then a rearrangement 
of the sum arising from (I5.7P being inserted into (15.61) yields formally 

oo 

— o = e r t - <«■<» 

1=1 trees T 
deg T=Z 

where the asterisk reminds us that the rules (Rl) and (R2) are being respected. 

According to the analysis above, the particular graphs appearing in the definition of 
deg T are the only possible single-node subgraphs of T proportional to a positive power 
of e. Since each tree is an analytic function of e, writing again ( • ) k for the fcth coefficient 
of the power series in e, we have 

oo oo 
rp _ £ fc rpk _ £ degT £ fc rpk+deg T 

fc=degT fc=0 

Hence, only trees with degree at most equal to I can contribute to b~iX l : 

^ £ =E T Ti = { T,* T ) ( 5 - 10 ) 



1=1 degT=l 



or, alternatively, 



S 2 X= ^*T + 0(e m ) (e->0) (5.11) 

deg T<i 

for each and every £=1,2,.... The expansion in (I5.9P is in fact just a compact way of 
writing (15. lip . We emphasize that the latter can be derived completely rigorously, for 
each value of £ separately, but resorting to the use of formal series allowed us to treat 
all orders of 82X at once. We call the series (I5.9P an asymptotic expansion of 82X] the 
partial sums ^de g r<£^ need not converge to 82X for any fixed e as t — > 00, but for 
a fixed I the error is bounded by an ^-dependent constant times |e| £+1 on the mutual 
domain of analyticity, |e| < e . 
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Example 5.1. The beginning of the asymptotic expansion ( 15. lip reads 

S 2 X = + 0(e 2 ) = + © + . + 



+ C(e 3 



5.2. Analyticity domain of trees. As already pointed out, all trees T above are 
analytic functions of (z, 9, e) on IT r x {|e| < e }. Due to the projections 5 2 appearing in 
( 15. 4p . they also satisfy T\ z=0 = d z T\ z=0 = 0, i.e., are elements of the space A± defined 
in ( 14. 7p . On this space, the inverse of /C = (q^ 2 ) ( see (14. 21) ) constructed in the proof 
of Lemma 14.31 satisfies 

r° ~ 

)C- 1 h(z,9)= K{s;z)h(ze 1 \9 + us)ds. (5.12) 

J — oo 

Consequently, we will now show that the analyticity domain of a tree in the z-variable 
is in fact much larger than the neighbourhood of [—1, 1] that is included in I1 T ; namely 
it includes the wedgelike region 

U T ,,? := {|z| < r} |J{arg2 G [-#,#] U [tt - tf, vr + tf] } C C 

(with a new r and "small" 

Lemma 5.2 (Analytic continuation of trees). Without affecting the analyticity domain 
with respect to e, there exist numbers < r < 1, < $ < tt/2, and < o < r\ such 
that each tree in the sums ( 15. 9p and ( 15.111) extends to an analytic function of (z, 9) on 
U T ,tf x {\%m9\ < a}. 

Proof. Observe that, as a polynomial, X<i is an entire function of z. On the other 
hand, $°(^) = 4arctanz = 2i(log(l — iz) — log(l + iz)), implying that |Qtn$°(;z)| < r\ 
in U T) ^ with r and •& sufficiently small. By Remark II. 3[ f(§°(z),9) is analytic, making 
the maps hS k ^ and X<i analytic on U T) ^ x {|$5tn^| < cr} for some < a < r/, where r] is 
determined by / and a by X< x (ultimately by / and uj). 

Suppose h = 82k is a map analytic on U Tj ^ x {|Qtn#| < a}. Then the integrand in 
(I5.12p is analytic in a neighbourhood of the latter set. By virtue of Fubini's theorem, 

(f lC- 1 h((,9)d(= [ I K(s;()h((e~< s ,9 + us)d(ds = 
Jr J -00 Jr 

for any closed contour T inside a sufficiently small neighbourhood of U T .^ and enclosing 
z. Hence, Morera's theorem yields analyticity of K~ l h with respect to z. As always, an- 
alyticity with respect to 9 follows from an exponentially decaying bound on the Fourier 
coefficients. Applying this argument at each node of a tree proves the claim. □ 

Proof of Theorem Since the number of terms in the sum in (15.101) is finite and the 
functions X<i and X° in X = X° + X<i + 5 2 X are analytic on U Tj ^ x {|Sm^| < a}, 
the analyticity of X follows from Lemma 15. 2[ 
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From the equations of motion, ( 11. lip , a Taylor expansion yields 

C 2 X = -C 2 X° + Q(X°) + V -±-D m n(X°)(Xf m , 

/ — ' ml 

m=l 

where the trigonometric degree of D m Q(X°) is N for e 7^ but vanishes at e = because 
X° does not depend on 6. For each k > 1, let n.^ stand for the trigonometric degree 
of Equating like powers of e in the expansion above, we infer two things. First, 
ri\ = N. Second, we must have, for each I > 2, 




n kx H h n km where k\ H h fc m = 4 

-/V + + h n fcm where ki H h & m = I — 1, 



because the trigonometric degree of a product is at most the sum of the trigonometric 
degrees of the factors; e iq (> e iq (> = e i2q9 and e^e"* 9 ' 9 = 1- 

Next, assume that n k < kN holds for each 1 < k < I — 1, recalling that this is the 
case if k — 1. Subsequently, the estimate for ri£ above becomes rig < £N. □ 
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